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SECTION  1 •  INTRODUCTION 


1.1  General 

NCAP  is  an  acronym  for  the  Rome  Air  Development  Center's 
Nonlinear  Circuit  Analysis  Program.  It  is  a  user-oriented 
computer  code  for  determining  the  nonlinear  transfer  functions  of 
electronic  circuits.  Py  utilizing  a  standard  set  of  circuit 
elements,  NCAP  can  analyze  networks  made  up  of  interconnections 
of  these  elements. 

NCAP  is  written  in  ANSI  Standard  FORTRAN  and  can  analyze 
networks  containing  up  to  approximately  500  nodes.  Sparse  matrix 
routines  are  used  to  decrease  core  storage  requirements  and 
increase  computational  efficiency  of  the  program. 

Structurally,  NCAP  solves  the  nonlinear  network  problem  by 
forming  both  the  nodal  admittance  matrix  (Y  matrix)  for  the 
entire  network,  and  the  first-order  generator  (current-source) 
excitation  vector,  for  all  of  the  linear  sources  in  the  entire 
network.  The  generators  can  be  located  at  any  node  in  the 
network,  and  can  have  any  desired  frequency,  amplitude  and  phase. 
Using  Gaussian  elimination  with  the  admittance  matrix  and  the 
current  vector  results  in  the  first-order  nodal  voltage  vector 
for  the  network,  the  elements  of  which  are  the  first-order 
transfer  functions  at  all  nodes  in  the  network  at  the  given 
excitation  freouency.  When  there  is  more  than  one  generator  at  a 
aiven  freouency,  the  first-order  transfer  function  will  be  the 
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total  transfer  function  due  to  the  superposition  of  the 
generators  since  the  first-order  transfer  function  is  a  linear 
function.  The  higher-order  transfer  functions  are  solved 
iteratively,  using  techniques  described  in  Volume  I,  Section  II. 

NCAP  is  relatively  simple  to  employ.  The  user  enters  a 
description  of  the  circuit  to  be  analyzed.  NCAP  interprets  the 
input  statements,  performs  the  nonlinear  analysis,  and  outputs 
the  results  in  printed  form.  There  are  several  types  of  input 
statements  which  are  needed  to  analyze  a  given  circuit.  The 
input  statements  define  the  topology  of  the  circuit,  the  circuit 
element  values,  the  linear  and  nonlinear  devices  used  in  the 
circuit,  the  circuit  excitation  and  the  order  of  the  analysis, 
the  desired  output,  the  data  modification  and  sweeping 
descriptions. 
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1.2  CPCANIZATIOK  OF  TUT  DOCUMENT 

This  document  consists  of  three  volumes.  Volume  I,  the 
Fncnneenng  Manual,  contains  introductory  remarks  on  the  history 
of  NCAP ,  a  review  of  the  Volterra  analysis  technique,  a 
Presentation  of  the  circuit  equations  programmed  in  NCAP,  ar.d 
detailed  discussions  of  the  active  NCAP  circuit  element  models  as 
well  as.  experimental  techniques  for  determining  the  model 
parameters.  Volume  II,  the  User's  Manual,  contains  an 
operational  introduction  to  the  NCAP  program  and  Volume  III,  the 
Programmer's  Manual,  .contains  narrative  descriptions  of  the  NCAP 
subroutines,  functional  flow  diagrams  depicting  the  general  flow 
of  loaic  through  the  program,  and  trace  maps  which  identify  the 
subroutine  linkages. 


1.3  PACKGROUND 


The  Rome  Air  Development  Center  (PADC)  has  been  tasked  with 
the  development  of  analysis  technioues  for  weapon  system 
electromagnetic  compatibility  (EMC )  assurance.  This  development 
effort  is  formally  known  as  the  Air  Force  Intrasystem  Analysis 
Program  (IAP) .  The  IAP  concept,  initiated  within  the  Air  Force 
in  the  early  1970's,  presently  consists  of  a  collection  of 
computer-aided  analysis  routines  for  addressing  various  aspects 
of  intrasystem  EMC.  The  analysis  routines  range  in  applicability 
from  an  overall  system  (aircraft,  satellite,  etc.)  level  EMC 
model  to  detailed  wire  coupling  and  circuit  analysis  models 
(NCAP) .  The  intent  is  to  provide  the  EMC  analyst  and  system 
designer  with  an  assemblage  of  technioues  for  solving  EMC 
problems  which  often  present  themselves  in  varying  degrees  of 
complexity.  NCAP,  under  the  general  framework  of  the  IAP,  could 
be  applied  after  a  more  coarse  analysis  has  indicated  a  potential 
EMC  problem  at  the  circuit  level.  At  this  point  NCAP  could  be 
used  to  examine  in  more  detail,  the  nonlinear  effects  which  are 
often  encountered  in  practice  and  can  severely  degrade  system 
performance.  It  can  be  shown  that  many  important  nonlinear 
interference  effects  such  as  desensitization,  intermodulation, 
and  cross  modulation  can  be  specified  in  terms  of  the  nonlinear 
transfer  functions,  functions  computed  by  NCAP.  A  knowledge  of 
the  magnitude  of  these  effects  could  be  very  valuable  in  the 
early  stages  of  system  design. 
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1.4  APPROPRIATE_ CLASS  QP_ PROBLEMS 

NCAP  is  applicable  to  a  specific  class  of  nonlinear 
problems.  Although  considerable  progress  has  been  made  in  the 
development  of  analytical  techniques,  there  is  no  general 
solution  technique  available  for  all  classes  of  nonlinear 
problems.  NCAP  is  devoted  to  the  class  of  nonlinear  problems 
known  as  "weakly  nonlinear".  The  nonlinearit ies  in  this  type  of 
problem  are  characterized  by  gradual,  rather  than  abrupt, 
variations.  Weakly  nonlinear  systems  encompass  an  extensive 
class  of  practical  systems.  They  include  small-signal  mixers  and 
square-law  detectors  as  well  as  the  large  class  of  "quasi-l inear" 
circuits.  By  quasi-l inear ,  we  mean  those  circuits  designed  to  be 
linear,  but  in  fact,  exhibit  nonlinear  behavior  over  some  region 
of  operation.  An  example  would  be  the  amplification  region  of  a 
transistor  amplifier.  Although  designed  to  linearly  amplify 
input  signals,  close  inspection  will  show  that  this  is  not  the 
case*  and  indeed,  the  output  wave  form  is  not  an  exact  replica 
of  the  input.  The  extent  of  this  "degradation"  is  often  of 
serious  concern  and  can  limit  amplifier,  and  consequently,  system 
performance. 
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THE  NONLINEAR  TRANSFER  FUNCTION  APPROACH 


SECTION  2^  THE_NONLINEAR_TRANSFER_  FUNCTION,  APPROACH 
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2.1  Introduction 

The  analysis  of  nonlinear  systems  is  complicated  by  the  fact 
that  no  single  analytical  tool  is  generally  applicable.  However, 
various  approaches  are  available  for  specific  classes  of 
nonlinear  systems.  For  the  special  case  of  weakly  nonlinear 
systems  the  Volterra  functional  series,  otherwise  known  as  the 
nonlinear  transfer  function  approach,  has  proven  to  be  an 
effective  analysis  method. 

Consider  a  weakly  nonlinear  circuit  with  input  x(t)  and 
output  y(t).  The  nonlinear  transfer  function  approach  models  the 
circuit  as  shown  in  Fig.  2.1.  This  model  consists  of  the 
parallel  combination  of  N  blocks  with  each  block  having,  as  a 
common  input,  the  circuit  excitation  x(t).  The  output  of  the  nth 
block  is  denoted  by  yn  ( t) ;  n»l,2,...,N.  The  total  system 

response  is  obtained  by  summing  the  outputs  of  the  individual 
blocks  so  as  to  yield 

N 

y(t)  =  y  x  ( t)  +  y2(t)  +  ...  +  yN(t>  -  z  yn(t).  (2*1“D 


The  nth  block,  characterized  by  the  nth-order  nonlinear  transfer 
function  Hn(fp  f2,  . ..,  f n )  ,  is  of  nth  order  in  the  sense  that 
multiplication  of  the  input  x(t)  by  a  constant  A  results  in 
multiplication  of  the  output  yn(t)  by  the  constant  An. 

Thus,  the  nonlinear  transfer  function  approach  represents 
the  output  of  a  weakly  nonlinear  system  as  a  sum  of  N  individual 
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responses.  The  linear  portion  of  the  circuit,  characterized  by 
the  conventional  linear  transfer  function  H^tf^),  generates  the 
first-order  component  of  the  response.  The  quadratic  portion  of 
the  circuit,  characterized  by  H  2  (f  1  generates  the 
second-order  component  of  the  response.  Additional  responses  are 

generated  in  a  similar  manner.  Blocks  above  Nth  order  arc  not 
included  in  the  model  because  it  is  assumed  they  contribute 

negligibly  to  the  output. 

A  weakly  nonlinear  system  is  completely  characterized  by  its 
nonlinear  transfer  functions.  Once  these  are  known,  the  system 
output  can  be  determined  for  any  given  input.  This  is  analogous 
to  the  case  of  linear  systems  where  knowledge  of  the  linear 
transfer  function  is  sufficient  to  completely  specify  the 
input-output  relation. 
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2.2  Sinusoidal  Steady-State  Response  of  Weakly  Nonlinear  Systems 

The  sinusoidal  steady-state  response  of  weakly  nonlinear 
systems  is  considered  in  this  section.  Let  the  input  in  Fig.  2.1 
consist  of  the  0  sinusoidal  tones  given  by 
Q 

x ( t)  =  Z  I  Xq  |  cos  (2nfqt  +  0q)  (2.2-1) 

q*i 

where  the  freouency,  amplitude,  and  phase  of  the  gth  tone  are 
denoted  by  f q ,  jxql,  and  e  ,  respectively.  In  order  to  express 
(2.2-1)  as  a  sum  of  complex  exponentials,  the  complex  amplitude 
of  the  qth  tone  is  defined  to  be 

Xq  -  I  X,  I  exp  tjeql.  (2.2-2) 


The  input  can  then  be  expressed  as 


x(t) 


1 

2 


Q 

E 

q=-Q 


Xq  exp( j2irfqt) 


(2.2-3) 


where  it  is  understood  that 


-q 


-f. 


-q 


=  x. 


=  0. 


(2.2-4) 


In  (2.2-4)  the  asterisk  denotes  complex  conjugate.  Whereas 
(2.2-1)  involves  only  positive  input  frequencies,  (2.2-3) 
suggests  an  interpretation  whereby  both  positive  and  negative 
frequencies  are  contained  in  the  input. 


r 


"tzzst 


1-13 


It  is  well  known  that  the  response  of  a  system  to  a  sum  of 
sinusoids  is  a  sum  of  sinusoids  whether  the  system  be  linear  or 
nonlinear.  However,  whereas  the  response  frequencies  of  a  linear 
system  are  identical  to  the  input  frecuencies,  new  frequencies 
not  present  in  the  inrut  are  generated  in  the  response  of  a 
nonlinear  system.  With  reference  to  (2.2-3)  ,  the  frequency 
content  of  x(t)  consists  of  the  2C  freouencies  • • •  '"f^ 
f  1 , . . .  ,f  f  q  .  The  output  of  the  weakly  nonlinear  circuit  of 
Fig.  2.1  contains  frequencies  of  the  form 


m_Q^_Q  +  +  rni^i  +  mQ^Q 


( mj-m_ f  ^  +•»•+ (mQ-m_Q ) f q 


(2.2-5) 


where  the  coefficients  •  k»-C , . . .  ,-l .  1 . . . .  #C  are  nonnegative 
inteaers.  For  a  specific  output  frequency  contained  in  the  nth  - 
order  portion  of  the  response  yn(t),  it  is  reouired  that 


+  m_^  +  + 


+  ■  n. 


(2.2-6) 


In  general,  the  freouencies  of  y  (t)  are  generated  by  all 
possible  choices  of  the  nonnegative  coefficients  m^  such  that  the 
constraint  given  by  (2.2-6)  is  satisfied. 

Fach  separate  choice  of  the  coefficients  r ^  in  (2.2-5) 
results  in  a  different  frecuency  mix.  To  characterize  the 
various  frecuency  mixes,  it  is  convenient  to  introduce  the 
frecuency  mix  vector 
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(2.2-7) 


m 


{m 


-Q 


•  •  •  l  I^q) 


For  example,  with  two  sinusoidal  inputs,  Q«2  and  m*(m_2»m_]  , 
m  j ,m  2  ) •  The  vectors  (1,0, 1,0)  (1,1, 1,0),  and  (0,0, 0,2)  then 
correspond  to  the  frequency  mixes  —  f l  l  ^  »  and 
( 2fj ) ,  respectively. 

when  the  input  x(t)  consists  of  0  sinusoidal  components,  the 
number  of  distinct  frequency  mix  vectors  m  associated  with  yn(t) 
is  given  by 


M 


( 2Q  +  n  -  1)  I 
~n1  r2Q~T~TTT 


(2.2-8) 


For  example,  with  0“2  (i.e.,  two  sinusoidal  inputs),  there  are 
four  distinct  frequency  mixes  associated  with  y 1 ( t ) ,  ten  distinct 
freauency  mixes  associated  with  y2(t),  and  twenty  distinct 
frequency  mixes  associated  with  y3(t).  The  various  frequency  mix 
vectors  and  the  corresponding  frequency  mixes  are  tabulated  in 
Table  1.  Note  that  the  elements  of  each  vector  sum  to  the  order 
n.  This  sum  is  referred  to  as  the  order  of  the  frequency  mix. 
For  example,  (2, 1,0, 2}  corresponds  to  a  frequency  mix  of  order  5. 

Frequencies  that  result  from  a  frequency  mix  are  termed 
intermodulation  frequencies.  The  same  intermodulation  frequency 
can  be  generated  by  several  different  frequency  mixes.  For 
example,  the  different  frequency  mixes  represented  by  (0,1, 0,2), 
(1,1, 0,3),  and  (0,2, 1,2)  all  generate  the  same  intermodulation 
freauency  2f2-f^.  For  clarity  of  discussion,  frequency  mixes, 
such  as  (-f  2  i  +3f  2  ^ '  are  enc^08ert  *n  parentheses  whereas 
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interirodulation  frequencies,  such  as  2f2-f^,  are  not. 

TABLE  2.2-1 

Pealization  of  {  n*_2 ,ir_1  ,ir2 )  for  n»l,2,  and  3 


*  2  » ^  ,r^  ,ir2 } 

Freouency  {n'_2,m_1<ir1,m2) 

Freouency 

Mix 

Mix 

n«l 

{1,0,0,01 

(-f2) 

{0,0,1,01 

<fl 

{0,1  .0,0} 

<-fl) 

{0,0,0,11 

<f2 

n»2 

{1,1, 0,0} 

(“«2  "fi> 

{0,0.1,11 

(fl+f2 

{1,0.1,01 

(-f2+fi) 

{0,1,0,11 

(-f  +f 

1  2 

{1,0,0.11 

(-f2  +f2> 

{0,1,1,01 

(-fi+fl 

{2,0,n,G} 

<-2f2) 

{0,0,0,21 

( 2  f  2 

{0,2  ,C  ,0} 

<-2f x) 

{0,0,2,01 

(2fi 

n»3 

{1,1,1.01 

{-f  2“f  j^  +  f  ) 

{0,1,1,11 

(”f l+f l+f 2 

{1,1.0,11 

(-f2~fl+f 2) 

{1,0,1,11 

("f2+fl+f2 

{2,1,0,01 

<-?f2-fi> 

{0,0,1,21 

( f  l+2f  2 

{0,2.1.01 

(“2f1+f1) 

{0,1,2,01 

<-fl*«l 

{1.2,0,01 

(-f2-2f1) 

{0,0,2,11 

(2fx+f2 

{2,0,0,11 

(”2f 2+f2) 

{1,0,0.21 

(”f 2+2f 2 

{2,0,1,01 

(-2f2+f1) 

{0,1,0,21 

<-fl+2f2 

{0,2,0,11 

(-2f1+f2) 

{1,0,2.01 

(-f2+2f1 

{3,0,0,01 

(-3f2) 

{0,0,0,31 

{ 3f  2 

{0,3,0,01 

<-3f1) 
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{0,0,3,01 

( 3  f  l 

From  (2.2-8)  it  is  obvious  that  M ,  the  number  of  distinct 
frequency  mix  vectors  associated  with  yn(t),  becomes  large  even 


for  moderate 

values 

of 

C  and  n . 

Fortunately, 

in 

most 

applications 

it  is 

not 

necessary 

to  evaluate  all 

of 

the  M 

resulting  responses.  In  general,  we  are  interested  in  only  those 
frequency  mixes  which  generate  intermodulation  frequencies  that 
fall  close  to  or  within  a  system's  pass  band.  Components  that 
fall  out  of  band  are  strongly  attenuated  and  are  not  usually 
cause  for  concern.  For  example,  consider  a  system  tuned  to  50 
MHz  with  a  1  MHz  bandwidth.  Assume  the  input  to  consist  of  two 
out-of-band  sinusoidal  tones  at  f  ^  *46MHz  and  f2ai48MHz.  The 
frequency  mix  vectors  associated  with  y^(t),  y2(t),  and  y^(t) 
are  listed  in  Table  1.  For  the  specified  values  of  f^  and 
f2  intermodulation  frequencies  in  y^(t),  arranged  in  numerically 
increasing  order,  are  -48,  -46,  46,  and  48  MHz.  All  of  these  are 
out  of  band  and  need  not  be  considered.  Similarly,  the 
intermodulation  frequencies  in  y2(t)  are  -96,  -94,  -92,  -2,  0,  2, 
92,  94  and  96  MHz.  These  also  fall  outside  of  the  system  pass 
band  and  need  not  be  considered.  Finally,  the  intermodulation 
frequencies  in  y3(t)  are  -144,  -142,  -140,  -138,  -50,  -48,  -46, 
-44,  44,  46,  48,  50,  138,  140,  142,  and  144  MHz.  Of  these  16 
intermodulation  frequencies  only  those  at  -50  and  50  MHz  fall  in 
band.  Hence,,  it  is  necessary  to  consider  only  the  two  responses 
of  y  3(t)  associated  with  the  frequency  mix  vectors  {2, 0,1,0}  and 
(0,1, 0,2}. 

The  next  step  is  to  evaluate  the  responses  corresponding  to 
the  frequency  mixes  of  interest.  The  nth-order  transfer  function 


HR(f i» . . . ,f n  )  plays  a  significant  role  in  determining  the 
amplitudes  and  phases  of  the  sinusoidal  components  in  the 


nth-ordei  portion  of  the  response  yn(t).  In  particular,  an 
nth-order  frequency  mix  characterized  by  the  frequency  mix  vector 
m  generates  the  response 

yn(t;5!>  '  lyJ  co»l2*f»t  +  V  •  -(2.2-9) 

From  (2.2-5)  the  response  frequency  is  the  intermodulation 
frequency 


i. 


m 


(ni1  -  m_1)f1  +  ...  +  (mQ  -  m-Q)fQ 


(2.2-10) 


The  response  amplitude  depends  upon  the  input  amplitudes  as  well 
as  the  magnitude  of  the  appropriate  nonlinear  transfer  function. 
Specifically, 

*  Ctn  (m_Q ! )  . . .  (m^lj  771  (mQ0 


+  m_1) 


m.Q> 


(2.2-11) 


•  •  • 


#  •  •  •  9 


where 

21-n  ,  f  /  0 

c  =  -  — 

m 

,  2"n  ,  f  =  0 

m 


(2.2-12) 

J 
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Observe  that  the  value  of  the  constant  C  depends  upon  whether  or 
not  the  interirodulation  frequency  is  zero.  The  exponents 
associated  with  the  input  amplitudes  are  specified  by  the  entries 
of  the  frequency  mix  vector.  The  freauency  mix  vector  also 
determines  the  n  arguments  of  the  nonlinear  transfer  function. 
Specifically,  the  frequency  fk  is  used  as  an  argument  m^  times* 
k—0, . . .  »-l  ,1  , . . .  ,0.  The  total  number  of  arguments  is  n  because 
of  the  constraint  in  (2.2-6).  The  order  in  which  the  n  arguments 
appear  is  actually  unimportant  because  the  nonlinear  transfer 
functions  are  symmetrical  in  their  arguments.  For  example, 


H3 (f i» f j, f 3)  =  H3(flff3,f2) 

=  H3(f2,f3,f1) 

-  H3(f3,f1,f2) 

*  H3(f3,f2,f1) 

(2.2-13) 

Finally,  the  response  phase  angle  is  a  function  of  the  input 
phase  angles  plus  the  phase  angle  of  the  nonlinear  transfer 
function.  In  particular, 

6  =  6  +  iL 

ym  m  vm  (2.2-14) 


where 


*ml  ”  ra-i)ei  +  +  (mg  -  m_Q)0Q, 


(2.2-15) 
. • ,fQ) 


Note  that  the  manner  in  which  the  input  phase  angles  combine  in 
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(2.2-15)  it  identical  to  the  way  in  which  the  input  frequencies 
combine  in  (2.2-10).  Also,  the  nth-order  transfer  function  in 
(2.2-15)  is  the  same  as  that  in  (2.2-11). 

For  a  specific  frequency  mix  vector,  yn(t;m)  is  readily 
evaluated  by  straightforward  application  of  (2.2-9)  through 
(2.2-15).  For  example,  consider  the  second-order  frequency  mixes 
characterized  by  the  vectors  (0,0, 0,2),  (0,0, 1,1),  and  (1,0, 0,1). 
The  corresponding  second-order  responses  are 

ya(t»(0,0,0,2))  -  |  |x2|2  |H2(trra)|  co»(2ir(2fa)t  ♦  2»a  ♦  *{o  a)  ] 

y2(t»{0, 0,1,1))  -  |X1 I  |x2|  |Ha(r1>Ca)|  co»(2w(f14f2)t  ♦  ♦  «2  ♦ 

*{0,0, 1,1)* 

(2.2-16) 

y2  (t;  (1, 0,0,’l) )  -  j  |X2|2  |H2(f2,-f2)|  cos  N>{lf0f0,i)  1 

In  general,  nonlinear  transfer  functions  associated  with  zero 
intermodulation  frequencies  are  real  quantities.  Hence, 

*{1  0  0  l)*s  either  zer0  degrees  180  degrees. 

It  is  important  to  remember  that  the  same  intermodulation 
frequency  can  be  generated  by  several  different  frequency  mixes. 
Consequently,  it  is  necessary  to  sum  the  responses  of  all  the 
frequency  mixes  which  generate  a  particular  intermodulation 
frequency.  The  total  response  at  frequency  f  is  denoted  by 

y(tjf).  To  illustrate  this  point,  consider  a  weakly  nonlinear 
circuit  for  which  the  highest  order  response'  is  of  fifth  order 
(i.e.,  N»5  in  Fig.  2.1),  Let  the  input  consist  of  two  sinusoids 

at  frequencies  fj^  and  f2  (i.e.,  C»2) .  The  output  contains  an 

intermodulation  frequency  at  2f  2  which  is  generated  by  the 

three  frequency  mixes  characterized  by  the  frequency  mix  vectors 


o 


(0.1,0,21,  (1,1,0,31.  ana  (0,2,1,21 

at  2fj-f^  is  given  by 


Hence,  the  total  response 


y(t)2f2  -  £j)  "  y3(t) (0,1,0,21)  +  y5 (t; (1, 1,0 , 3l)  +  y5 (t; (0, 2, 1,2) ) 

*  T  l*l!  I X2  I  1 0-j ^2 * ^2^  I 

coa(2n(2f2-£1)t  +  28j  -  0X  +  *{0(lj0,2)! 

+  |xxl  |x2|  |  Hg  (-f  2 , -f  2 ,  f  2»  f  2 ,  f  2)  | 

cos  [2ir  (2f  2  -  f^t  +  202  -  0j.+  ♦{!, 1(0 ,3) 1 

+  — y  | X^ | ^  |X2 I2  I Hg (-f 3 ,-f  3,f 3 , f 2' £2^ I 

cos  (2ir  <  2£  2  -  +  2«2  -  01  +  *(0,2,1,2)1 

(2.2-17) 


O 


In  general,  the  angles  *(0_ li0_2] ,  *{1, 1(0, 3)  ,  and  1{o,2,i,2)  «e 
not  equal.  The  various  terms  in  (2.2-17)  combine  as  shown  in  the 
phase  diagram  of  Fig.  2.2.  Note  that  both  magnitude  and  phase 
angles  of  the  nonlinear  transfer  functions  are  important  in 
determining  the  total  response  at  a  particular  intermodulation 
frequency.  For  very  small  values  of  lxj.1  and  |x2l,  the 
response  is  dominated  by  the  third-order  term  and  the  fifth-order 
terms  may  be  ignored.  However,  the  fifth-order  terms  become 
significant  as  the  input  tones  are  increased  in  amplitude. 
Finally,  for  very  large  values  of  J  x^^  |  and  .J  X2  I ',  the  response  is 
dominated  by  the  fifth-order  terms. 
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Having  demonstrated  the  role  played  by  the  nonlinear 
transfer  functions  in  the  sinusoidal  steady-state  response  of 
weakly  nonlinear  circuits,  the  procedure  used  by  NCAP  to 
determine  the  nonlinear  transfer  functions  is  discussed  next. 


y 


\ 


Fig.  2.2  Phasor  diagram  illustrating  total  response  at  the 
intermodulation  frequency  2f2-fi. 
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2*3  NCAP  Procedure  for  Determination  of  H  (f  ^  f  3, . . .  ,f  ) 

The  basic  nonlinear  circuit  elements  founa  in  NCAP  are  the 
nonlinear  resistor,  nonlinear  capacitor,  nonlinear  inductor,  and 
nonlinear  controlled  sources.  These  are  adequate  for  modeling 
the  nonlinearities  of  most  electronic  devices  such  as  diodes, 
vacuum  tubes,  transistors,  integrated  circuits,  etc.  In  weakly 
nonlinear  circuits  electronic  devices  are  typically  operated  over 
a  localised  region  of  their  characteristics.  Therefore,  it  is 
possible  to  expand  the  nonlinearities  in  power  series  about  the 
quiescent  operating  point  in  terms  of  incremental  variables.  The 
series  representations  of  the  basic  nonlinear  circuit  elements 
are  presented  next. 

Let  er(t)  and  ir  (t)  denote  the  incremental  voltage  and 
current  associated  with  a  nonlinear  resistor.  The  incremental 
current  through  the  resistor  may  be  expressed  as 


i  (t)  -  E  g.(er(t)]k 

r  k-1  K  r 


(2.3-1) 


Note  that  g^  is  the  conductance  of  the  linear  resistor  normally 

included  in  a  linear  incremental  equivalent  circuit,  similarly, 

let  e  (t)  and  i  (t)  denote  the  incremental  voltage  and  current 
c  c 

associated  with  a  nonlinear  capacitor.  The  incremental  current 
through  the  capacitor  is  given  by 


ic<t) 


CKtao(t»> 


deo(tl 


(2.3-2) 


) 
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The  term  corresponding  to  k«0  is  in  the  usual  form  for  a  linear 
capacitor.  Hence,  CQ  is  the  capacitance  of  the  linear  capacitor 
to  be  inserted  into  a  linear  incremental  equivalent  circuit.  '*  For 
a  nonlinear  inductor  the  series  representation  of  the  current  is 

i « ( t)  -  £  rkt;^  e,  (zf  dzjk  (2.3-3) 

where  e  t  (t)  and  i£(t)  are  the  incremental  variables  associated 
with  the  inductor,  r ^  is  the  reciprocal  inductance  of  the  linear 
inductor  that  would  be  placed  in  a  linear  incremental  equivalent 
circuit. 

The  final  nonlinear  circuit  elements  to  be  discussed  are  the 
nonlinear  controlled  sources.  These  occur  frequently  in 
electronic  device  models  and  are  two-terminal  elements  whose 
terminal  voltage  or  current  is  a  nonlinear  function  of  either  the 
control  voltage  ev  or  the  control  current  i  from  some  other  part 
of  the  circuit.  If  the  control  variable  pertains  to  a  particular 
circuit  element  X,  this  element  is  called  the  controlling 
element.  The  series  representations  for  the  four  types  of 
nonlinear  controlled  sources  are: 


a)  Voltage  -  controlled  voltage  source 

ec.(t)  -  "k'V4”1'  - 

b)  Current  -  controlled  voltage  source 


(2.3-4) 
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(2.3-5) 


«cs(t)  -  R;1  rksUx(t)1 


c)  Voltage  -  conttolled  cuttent  soutce 


ics(t)  *  k“x  *ks"-x 


2  9uJe  (t))K  , 


(2.3-6) 


d)  Cuttent  -  controlled  cuttent  soutce 

lcs(t)  *  *  akUx(t)lk  *  (2.3-7) 

k  =  l 

The  linear  term  in  each  expansion  defines  the  lineal  conttolled 
soutce  to  be  used  in  linear  incremental  eouivalent  circuits.  The 
linear  coefficients  ,  r^,  gls,  and  die  referred  to  as  the 
voltage  amplification  factor,  the  mutual  resistance,  the  mutual 
conductance,  and  the  cuttent  amplification  factor,  respectively. 

The  procedure  employed  by  NCAP  to  solve  for  the  nonlinear 
transfer  functions  makes  use  of  nodal  analysis.  This  approach 
was  chosen,  as  opposed  to  loop  analysis,  because  electronic 
circuits  tend  to  have  many  fewer  nodes  than  loops.  Assume  the 
circuit  to  be  analyzed  contains  K+l  nodes.  One  of  the  nodes, 
usually  the  ground  connection,  is  assigned  to  be  the  reference  or 
datum  node.  Node-to-datum  voltages  ate  then  defined  for  each  of 
the  remaining  K  nodes.  Let  v  ^  (t)  denote  the  node-to-datum 

voltage  for  node  j*  j-l,2,...K.  (A  supet script  is  used  to  avoid 
confusion  with  subscripts  which  are  reserved  fct  the  designation 
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cxau  [ 


of  order.)  The  nth-order  transfer  function  at  the  jth  node  is 
denoted  by  H^(f . #fn) .  Hence,  the  jth  node  has  associated 
with  it  a  set  of  nonlinear  transfer  functions  of  order  1  to  N. 
The  excitation  is  assumed  to  be  a  voltage  source  vg(t)  with 
generator  impedance  zg(f).  For  convenience  of  discussion,  let 
the  source  be  connected  between  node  1  and  the  reference  node. 

The  first  step  in  the  procedure  is  to  write  Kirchhoff's 
current  law  at  each  of  the  K  nodes,  other  than  the  reference 
node.  The  current  through  each  nonlinear  circuit  element  is  then 
expressed  in  a  power  series  where  the  series  variables  are 
node-to-datum  voltages.  The  linear  term  in  each  series  can  be 
represented  by  a  linear  circuit  element.  These  are  lumped  with 
the  linear  portion  of  the  network  to  form  the  augmented  linear 
network.  In  terms  of  the  K  Kirchhoff  current  law  equations,  this 
latter  step  is  accomplished  by  rearranging  the  equations  such 
that  only  terms  linear  in  the  node-to-datum  voltages  appear  on 
the  left  hand  side.  The  right  side  of  the  equations  then  contain 
only  the  source  term  and  the  nonlinear  terms  from  the  power 
series  expansions  of  the  nonlinear  circuit  elements. 

The  solution  procedure  is  recursive  in  nature.  The 
augmented  linear  circuit  is  first  excited  by  allowing  vg(t)  to  be 
the  single  complex  exponential  exp[  j2rrf1t  ].  This  permits 
determination  of  the  linear  transfer  function.  Then  the  sum  of 
two  complex  exponentials,  given  by  Vg(t)  ■  exp[j2nf^t  ]  + 
exp[  j2nf2t  ],  is  applied.  This  yields  the  second-order  transfer 
function  in  terms  of  the  linear  transfer  function.  The  procedure 
continues  with  one  additional  complex  exponential  being  added  to 


the  input  et  each  step.  Therefore*  at  step  n ,  the  input  consists 
of  the  sum  Vg(t)»exp(  j2*fjt  1  ♦  ...  ♦  exp(j2»fnt  ).  The  nth  -order 
transfer  function  is  then  found  to  be  constructed  from  all  cf  the 
previously  determined  lower  order  nonlinear  transfer  functions. 

Let  the  admittance  matrix  of  the  augmented  linear  network  be 
denoted  by  (Y ( f ) ] •  Also*  define  the  nth-order  nonlinear  transfer 
function  vector  to  be 


V£1 . (n> 


Hn<fl . *n>l 


Hn(fi* ...,fn)  | 


Hn(fi*...*fn> 


(2.3-8) 


The  first-order  transfer  functions  for  the  various  node-to-datum 
voltages  are  simply  the  conventional  linear  transfer  functions 
and  are  obtained  by  solving  the  matrix  equation 


l  Y(ft)  ]  Hx (f x) 


r  1 

‘  w 


0  0 


(2.3-9) 


where  T  denotes  the  transpose  operator. 

The  second-order  nonlinear  transfer  functions  are  determined 
next  and  the  procedure  continues  recursively  for  each  of  the 
higher  transfer  functions.  In  general*  the  nth-order  transfer 
functions  are  solutions  to  the  matrix  equation 


O 


I  If ff ,  +  ...  +  f  )  J  H  (f. . f  ) 

x  n  — n  x  n  • 
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(2.3-10) 


where  the  nth-order  source  vector,  l_(f . , . . . ,f  ) ,  is  a  function 

-n  l  n 

of  the  lower  order  nonlinear  transfer  functions.  At  each  stage 
of  the  solution,  note  that  it  is  the  admittance  matrix  of  the 
augmented  linear  network  which  is  evaluated  at  the  appropriate 
frequency  and  then  inverted. 

The  jth  entry  in  the  nth-order  source  vector  arises  from  the 
nth-order  currents  associated  with  the  nonlinear  circuit  elements 
connected  to  node  j.  Evaluation  of  these  entries  becomes 
increasingly  tedious  as  n  grows  in  value.  The  situation  is  aided 
by  means  of  a  recursion  relationship  which  hr:  ’  .en  developed  for 
this  purpose.  The  exact  form  of  the  relationship  depends  upon 
whether  the  nonlinearity  is  a  zero-memory  nonlinearity  (i.e.,  a 
nonlinear  resistor  or  a  nonlinear  controlled  source),  a 
capacitive  nonlinearity,  or  an  inductive  nonlinearity.  The 
zero-memory  nonlinearity  is  discussed  first. 

Consider  a  zero-memory  nonlinearity  of  the  form 


i(t)  -  E  a.[vp{t)]k 

k-2  *  (2.3-11) 

where  v  p  (t)  denotes  the  node-to-datum  voltage  for  node  p.  The 
linear  term  has  been  omitted  in  (2.3-11)  because  it  is  included 
with  the  augmented  linear  network.  If  the  nonlinear  circuit 
element  is  connected  to  node  j,  the  jth  entry  of  the  nth-order 
source  vector  is  given  by 


(2.3-12) 


I-J  (fr..-,fn)  -  vn 

where  the  overbar  is  used  to  denote  the  arithmetic  average  of  the 
n!  terms  generated  by  all  possible  permutations  of  the  n 
frequencies.  For  example, 


X3  (,X',2'f3)  *  V3  ,fl'f2'f3>  -  f  tv3  (fj.fj.fj) 

+  *{  +  Vl  (t2,tytj)  +  (fj.fj.fj)  +  V|  (fj.fj.fj)] 


The  quantity  v  ■}  (fi»...,fn)  may  be  expressed  as 


(f1#...,fj  «  z 


*-2  ak  . f„> 


(2.3-13) 


(2.3-14) 


where  a  recursion  relationship  for  H  p  ,(f. )  is 

n#x  i  n 


n-k+1 


(2.3-15) 


Recall  that  the  ith-order  transfer  function  at  the  pth  node  is 
denoted  by  (f1,...,fi).  (f1#.;.,fn)  has  the  properties 
that 


H*j  (f •  #fy^) 


I 


Hntfl" 


.,fn)  ;  k«l 


}  k>n  •  . 


(2.3-16) 


1-29 


Because  of  the  permutations  implied  by  the  overbar  in  (2.3-12), 
the  order  in  which  the  n  frequencies  are  assigned  H  P  (  )  and 
H  p  .  ,  .  .  (  )  in  the  sum  of  (2.3-15)  is  immaterial  as  long  as 

each  frequency  appears  once  and  only  once.  For  the  sake  of 
consistency,  let  the  frequencies  •  be  assigned  consecutively 
beginning  with  the  factors  requiring  the  fewest  arguments. 

By  way  of  example,  the  following  results  are  obtained 
through  application  of  (2.3-15)  and  (2.3-16): 

«?,!«!>  *  "?<£1> 

H5,l(fl,f2>  "  H2(fl'f2)  '  “  Hl(fl)  Hl(f2) 


H3,l(£l'f2'£3>  *  , 


-  Hlt£l>  Hlt£2>  Hl(£3> 


H4,l(fl'f2'f3'f4}  "  H4(fl,f2,f3,f4)  '  "  ^H2<fl#f2>  ”2^3'^ 


+  2Hj(f1)  HP(f2,f3,f4)]  , 


H4,3(£l'£2'£3<£«>  *  »!<«!>  «?<£2>  «?«3-£4>  ' 

H4,4(fl'f2'f3'f4}  "  H?(fl)  Hl(f2)  Hl<f4> 

(2.3-17) 

From  (2.3-12)  and  (2.3-14),  it  follows  that  the  jth  entries  of 
the  nth-order  source  vector  for  n«2,3,  and  4  are 
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-  “2  ^  <£1>  «?  (f2> 

I^(firf2,f3)  -  2^  H^)  HP(f2,f3)  ♦  a3  ^(ty  HP(f2)  HP(f3) 

_  _  (2.3-18) 

i2(frf2,f3,f4>  -  +  2HP(fx)  nP(f2,f3,f4)] 

♦  3a3  ri[(f2)  nP(f3,f4)  +  a4  i^(f2)  nP(f3)  nP(f4) 

Observe  that  the  nth-order  entry  involves  only  nonlinear  transfer 
functions  of  lower  order. 

A  nonlinear  capacitor  is  now  assumed  to  be  connected  between 
nodes  j  and  p.  Let  the  nonlinear  portion  of  the  current  through 


the  capacitor  be  given  by 


I  c.  [v^t)] 

k-1  K 


k  dv^p(t)  d 

nr  *  % 


'  £  HT  [vjp(t))lt'  (2.3-j 

k«2 


where,  as  in  (2.3-11),  the  linear  term  has  been  omitted.  The 
voltage  across  the  capacitor  is 


v*(t)  *  v^(t)  -  vP(t)  . 


(2.3-20) 


The  jth  entry  in  the  nth-order  source  vector  is  given  by 
(2.3-12).  However,  now 

. fn>  -  j2»(f1+  ...  +  fnl  S  HjPk(f1 . fn>  (2.3-21) 


where 


Hn^k*fl ,mm"V  "  Hn,k^fl#"*,^n^  “  Hn,k^fl'***'fn^  *  (2.3-22) 
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The  recursion  relation  in  (2.3-15)  may,  once  again,  be  used  with 
the  superscript  p  replaced  by  jp.  For  example,  application  of 
(2.3-12),  (2.3-21),  and  (2.3-15)  for  n-3  results  in 

I3(fl'f2'f3)  “  j2w(fl  +  f2  ♦  f3}  [2  4  H1P  (fl>  H2P  {f2'f3} 

(2.3-23) 

♦  hJp  (fx)  hJp  (f2)  hJp  (f3)]  .* 

The  differences  between  the  expressions  in  (2.3-18)  and  (2.3-23) 
for  I  ^(f1,f2,f3)  should  be  noted  carefully. 

Finally,  assume  a  nonlinear  inductor  is  connected  between 
nodes  j  and  p.  The  nonlinear  portion  of  the  current  through  the 
inductor  is 

i(t>  -  J  fk  t'l  v*  (z>  dz]*  .  (2.3-24)  ( 


As  with  the  zero«*memory  nonlinearity  and  the  nonlinear  capacitor, 
the  jth  entry  in  the  nth-order  source  vector  is  given  by 
(2.3-121.  However,  now 


n 

E 

k-2 


G&(fl . £n> 


(^.3-25) 


jp 

where  a  recursive  relationship  for  G  n  k  1' • • • » 1S 


<&(fl . fn» 


n-k+1  hJp  (f1# _ ,f±) 

i«l  j27r(f1  +  ...  +  fA) 


4i,k-l,fi+lMM,fn) 


(2.3-26) 
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G  3P.  (£,,.. .,f  )  has  the  properties  that 
n,K  1  n 

/  H^P  (f1#...#fn) 


j2w(f1  fn)  ;  k»l 


G^P  (f  f  ) 


(2.3-27) 


;  k>n  . 


As  with  (2.3-15)  ,  let  the  n  frequencies  in  the  sur,  of  (2.3-26)  be 
assigned  consecutively  beginning  with  the  factors  recuirino  the 
fewest  arguments. 

As  an  illustration  of  (2.3-26)  and  (2.3-27),  note  that 


GM(fl) 


*1  (fj) 


fl(fwfj)  ,  gJp  2(frf2) 

2,1  1  2  j2n(f1  +  f2)  2'2  1  2 


HlP(fl)  HlP<f2) 
j2TTf1  j27Tf2 


Id  ^3^  ^l'^2'^ 

1  1  J  j27r(f1  +  f2  +  f3) 


,  G3  2(fl'f2'f3)  »  2 


j2irf. 


H^P(f2,fO 


,  GJP3(f1,f2,f3) 


1P<£1> 

Hip(f2> 

H^p(f3) 

j2*f1 

j2irf  2 

j27rf3 

j2ff(f2  +  f3)  ‘  J  j2irf1  j2irf2  j2irf3 

(2.3-28) 

The  jth  entry  in  the  nth-order  source  vector  is  obtained  by  a 
straightforward  application  of  (2.3-12)  and  (2.3-25).  For 
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example,  with  n«3 


2  j27Tf1  j2lT(f2  ♦  f3) 


♦  r. 


H^P(fx)  H|P(f2)  HjP(f3) 

j2irf1  j2ir£2  j2irf3 


(2.3-29) 


*e  see,  therefore,  that  the  entries  for  Iff, , . . . ,f  )  are  readily 
obtained.  The  key  step  in  the  procedure  involves  use  of  the 
recursion  relations  given  by  (2.3-15)  and  (2.3-26). 

The  procedure  used  by  NCAP  to  determine  the  nonlinear 
transfer  functions  is  now  illustrated  by  means  of  an  example. 
Consider  the  weakly  nonlinear  incremental  circuit  shown  in  Fig. 
2.3. 


Fig.  2.3  Weakly  nonlinear  incremental  circuit. 


Enclosed  in  boxes,  to  distinguish  them  from  linear  elements,  are 
a  nonlinear  capacitor,  a  nonlinear  resistor,  and  a  nonlinear 
inductor.  The  circuit  also  contains  a  nonlinear 
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Equivalent  circuit  of  nonlinear  network  shown  in  Fig. 


voltage-controlled  current  source  which  is  designated  by  iQff  (t)  . 
The  series  expansions  for  the  incremental  currents  through  the 
nonlinear  elements  are 


V*) 


CjJv^t)  -  v2(t)]k*  ^(t)  -  y2(t) 


ir(t)  -  2  gk[v2(t)]k 

k*l 

m 

it(t)  -  2  rk{^[v2(s)  -  v3(*)]ds)k 


(2.3-30) 


Wv1<«ik  • 

For  nodal  analysis  it  is  convenient  to  convert  the  independent 
voltage  source  vg  (t)  into  a  equivalent  current  source  i N< t)  by 
means  of  Norton's  theorem.  In  addition,  separating  each  power 
series  in  (2.3-30)  into  a  linear  term  plus  second-and 
higher-order  terms,  each  nonlinear  circuit  element  may  be 
replaced  by  a  linear  circuit  element  in  parallel  with  a  nonlinear 
voltage-controlled  current  source.  This  is  shown  in  Pig. 4  for 
the  circuit  of  Pig.  3.  With  this  interpretation,  the 
nonlinearities  manifest  themselves  as  controlled  sources  driving 
the  augmented  linear  network.  In  Pig.  4 

v12(t)  -  v1(t)  -  v2(t) 

(2.3-31) 

v23(t)  -  v2(t)  -  v3(t)  . 
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□ 


The  admittance  matrix  of  the  augmented  linear  network  is 


[Y(f)J 


-  +  g  +  j2iTfcn 

Zff)  a  o 


-j2irfCr 


-j27rfC0 
gx  +  j27rfC0  + 

.  JL 

j2irf 


j2irf 


ll 

j2irf 


j27TfJ 

(2.3-32) 


The  linear  transfer  functions  are  obtained  by  solving  the  matri> 
equation  (2.3-9)  while  the  higher  order  nonlinear  transfer 
functions  are  solutions  to  the  matrix  equation  (2.3-10).  For  n»2 
and  3  the  nth-order  source  vectors  are  given  by  (2.3-34)  where 


Hi2(fl . fn> 


aj(f fn)  -  Hj(f1M..rfn) 


Hn3*fl'”*'fr>' 


Hn(fl'“* ,fn)  ’  Hn(fl',,,,fn) 


(2.3-33) 


1-U7 

\ 

\ 


-j2ir(f^  ♦  f2)  ItJ2^)  "i2^^ 


-  {  -*2 


3  2  «?3«i>  »»i3(fj> 

W  W  ”r2  T*n^  15t?"  +  ^(fi  * 


C1 

f2}  ~r 


n^fr)  ) 


H?^(f.)  nf3{fj  -r - : — 

{  r2  1T3”  153J  *  *u  *1  (fl>  Vty  } 


(2.3-34) 


-j2ir(f1  ♦  t j  ♦  f3)  (2  ^  ♦  p-  H32^)  H“(f3)J 


{-[2g2  H2^)  1 

£j<  Wj>  . 

H?3(f.) 

’t2r2  -^fTf 

H^3(f2,f3) 

jtor(f2*r3) 

♦  32*«x  ♦  t2 

*  fJ>  «  7 

(l2r2  “73w?“ 

H^3(f2,f3) 

jwfjifj) 

♦I  ^2,  H^)  H3(«j.«3)  ♦  13  "iW  Hl<f3),) 
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Note  that  the  jth  entry  in  the  Nth-order  source  vector  consists 
of  the  sum  of  the  contributions  froir  the  currents  through  each  of 
the  nonlinear  circuit  elements  connected  to  node  j.  Also,  the 
sign  of  each  contribution  is  positive  or  negative  depending  upon 
whether  the  current  flows  into  or  out  of  the  jth  node. 
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2.4  Introduction  to  Electrical  Device  Modeling 

The  nonlinear  analysis  of  an  electronic  circuit  requires 
that  each  electronic  device  be  replaced  by  an  equivalent 
nonlinear  circuit  model.  Most  of  the  nonlinearities  encountered 
in  electronic  devices  are  capable  of  being  modeled  by  various 
combinations  of  nonlinear  resistors,  capacitors,  inductors, 
and/or  controlled  sources.  Several  different  circuit  models  are 
typically  available  for  each  device.  Some  are  quite  elaborate 
while  others  are  relatively  crude.  Models  which  are  highly 
accurate  over  large  ranges  of  amplitude  and  frequency  tend  to  be 
rather  complicated  and  computationally  inefficient-  Furthermore, 
their  complete  description  frequently  entails  a  considerable 
amount  of  effort  in  order  to  obtain  values  for  the  large  number 
of  model  parameters  typically  required.  Less  sophisticated 
models  tend  to  be  more  efficient  and  easier  to  specify.  However, 
since  they  arise  by  either  neglecting  or  approximating  various 
physical  effects,  they  are  limited  in  their  usefulness. 
Nevertheless,  simpler  models  can  provide  acceptable  results  in 
many  situations. 

NCAP  provides  stored  nonlinear  circuit  models  for  the 
semiconductor  diode,  the  bipolar  junction  transistor  (3JT) ,  the 
junction  field-effect  transistor  (JFET) ,  the  vacuum-tube  diode, 
the  vacuum-tube  triode,  and  the  vacuum-tube  pentode.  Although 
the  models  used  are  not  the  most  sophisticated,  they  have  been 
shown  to  yield  good  results  in  practice.  The  circuit  models  are 
presented  in  the  following  sections,  where  possible,  analytical 
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expressions  based  upon  the  physics  of  the  devices  are  given  for 
the  nonlinearities.  Otherwise,  empirical  relations  are  utilized. 
The  parameters  associated  with  each  device  are  defined.  Finally, 
power  series  expansions  of  the  incremental  currents  through  the 
nonlinear  circuit  elements  are  presented. 

The  operation  of  electronic  devices  is  usually  depicted  in 
terms  of  their  static  characteristics.  These  are  clots  relating 
a  device's  terminal  voltaoes  and  currents  as  obtained  from 
measurements  performed  under  dc  or  very  slowly  varyina 
conditions.  Therefore,  the  static  characteristics  of  actual 
devices  do  not  include  the  effects  of  energy-storage  elements  and 
frequency-dependent  parameters.  Hence,  the  approach  used  to 
develop  the  incremental  nonlinear  eauivalent  circuit  models  was 
to  first  obtain  a  dc  resistive  equivalent  circuit  to  model  the 
static  characteristics.  Capacitors  and  inductors  were  then  added 
at  one  or  more  strategic  locations  to  account  for  the  device's 
frequency  behavior. 

The  region  of  operation  for  an  electronic  device  is 
established  by  employing  a  dc  biasing  circuit  which  determines 
the  dc  or  quiescent  operating  point.  The  application  of 
additional  signals  then  results  in  the  device  voltages  and 
currents  varying  in  some  neighborhood  about  the  operating  point. 
It  is  important  to  distinguish  between  total,  dc,  and  incremental 
variables.  In  this  manual  total  variables  are  denoted  by  lower 
case  symbols  with  upper  case  subscripts.  For  example,  iD  and 
vD  represent  the  total  current  and  voltage,  respectively,  of  a 
diode.  The  dc  value  of  a  variable  is  denoted  by  an  upper  case 


symbol  with  an  upper  case  subscript.  For  a  diode  the  dc  current 
and  voltage  are  denoted  by  ID  and  V  D ,  respectively.  Incremental 
variables  are  represented  by  lower  case  symbols  with  lower  case 
subscripts  and  are  defined  to  be  the  difference  between  the  total 
and  dc  variables.  For  the  diode  the  incremental  current  and 
voltage  are  given  by 


(2.4-1) 


In  the  following  sections  both  global  circuit  models,  involving 
total  variables,  and  incremental  circuit  models,  involving 
incremental  variables,  are  presented.  In  the  schematic  diagrams 
boxes  are  places  around  resistors  and  capacitors  to  denote 
nonlinear  circuit  elements. 
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2.5  Semiconductor  Diode  Model 


A  global  model  for  the  semiconductor  diode  is  shown  in  Fig. 
2.5.  (D 


Fig.  2.5  Global  model  for  semiconductor  diode. 


The  four  circuit  elements  in  the  global  model  are: 

1)  R,  the  junction  leakage  resistance. 

2)  R j t  the  nonlinear  resistor  due  to  the  diode  junction. 
Its  current-voltage  relationship  is  given  by 


ij  ■  i,  t«p<  —  >  -  i] 

nW 


(2.5-1) 


3)  CBt  the  nonlinear  transition  capacitance  of  the 
depletion  layer ,  also  known  as  the  barrier  capacitance  or  the 
space-charge  capacitance.  It  is  given  by 


K  (-VD)~U 


9 
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(2.5-2) 


4)  Cd,  the  nonlinear  diffusion  capacitance  whose  value 
depends  linearly  on  the  junction  current  in  accordance  with  the 
relation. 

Cd  "  Cj  LJ  (2.5-3) 

The  parameters  which  appear  in  (2.5-1),  (2.5-2),  and  (2.5-3)  are 
defined  as  follows? 

Ig*  diode  saturation  current, 

-19 

q  •  electron  charge  -  1.6  x  10  coulombs, 

-23 

k  »  Boltzmann's  constant  ■  1.38  x  10  joules/degree  Kelvin, 
T  ■  junction  absolute  temperature  in  degrees  Kelvin, 

?.  n  *  diode  nonideality  factor, 

K  ■  transition  capacitance  at  1  volt  reverse  bias, 
v  m  junction  grading  constant, 

t 

cj»  diffusion  capacitance  constant. 

Near  room  temperature  (T»290#  K)  kT/q  is  approximately  25 
millivolts. 

The  global  model  may  be  simplified  when  the  diode  is  either 

forward-or  reverse-biased.  We  first  consider  forward-biased 

conditions  for  which  vk>0.  Since  the  forward  resistance  of  Rj  is 

then  much  smaller  than  that  of  R,  negligible  current  flows 

through  R  and  the  leakage  resistance  may  be  ignored.  Also,  when 

the  diode  is  forward  biased,  vD  is  essentially  constant  relative 

to  the  junction  current  ij  .  Therefore,  the  transition 

capacitance  C_  may  be  treated  as  a  constant  which  is  denoted  by 
Cj.  A  global  model  for  the  diode,  in  forward  bias,  is  shown  in 


c 
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f 


0 
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Fig.  2.6. 


Fig*  2.6  Forward-biased  global  model  Cor  the  semiconductor 
diode  (vD  >  0) . 


The  incremental  current  through  R ^  may  be  written  as 

m 


v°  -  & 


(2.5-4) 


where 


nkT 


(2.5-5) 


and 


%  "  JcIJJ  gl  gk-l 


(2.5-6) 


Recall  that  ID  is  the  dc  diode  current.  The  linear  incremental 
resistance  of  the  diode,  denoted  by  r^,  is  the  reciprocal  of  the 
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^  i 


conductance  g  L<  Note  that  it  inversely  proportional  to  ID  . 
The  series  expansion  for  the  incremental  current  through  C^is 


w« 


E 

lo— 0 


'•'dk 

Jc+T 


d 

at 


{  [v,,(t)]k+1  } 


(2.5-7) 


where 


(2.5-8) 


and 


(2.5-9) 


The  coefficients  in  the  expansion  for  the  diffusion  capacitance 
are  seen  to  be  intimately  related  to  those  in  the  expansion  for 
Rj.  The  forward-biased  nonlinear  incremental  equivalent  circuit 
for  the  semiconductor  diode  is  shown  in  Fig.  2.7.  Note  that  only 


Fig.  2.7  Forward-biased  nonlinear  incremental  equivalent 
circuit  for  semiconductor  diode. 

incremental  variables  appear  in  the  model.  Nevertheless  r^,  Cdo 
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and  the  aeries  coefficients  are  functions  of  the  quiescent 
operating  point.  The  linear  incremental  equivalent  circuit 
consists  of  rd,  Cj»  and  cdQ  all  in  parallel. 

The  diode  is  reverse  biased  when  the  diode  voltage  is 
negative.  The  junction  current  ij  and  the  diffusion  capacitance 
Cd  becomes  negligibly  small  when  vD  <0.  Therefore.  under 
reverse-biased  conditions,  the  diode  global  model  is  approximated 
by  the  circuit  shown  in  Pig.  2.8.  Cs  is  now  the  only  nonlinear 
element  contained  in  the  model . 


Pig.  2.8  Re  verse -biased  global  model  for  the  semiconductor 
diode  (>b<0) . 


(2.5-10) 


(2.5-11) 


and 


CaOc-l) 


(2.5-12) 


V  ^ ,  of  course,  is  the  dc  diode  voltage.  The  reverse-biased 
nonlinear  incremental  equivalent  circuit  for  the  semiconductor 
diode  is  shown  in  Fig.  2.9.  In  practice,  only  a  finite  number  of 
terms  from  the  infinite  sum  are  needed. 


Fig.  2.9  Reverse-biased  nonlinear  incremental  equivalent  circuit 
for  semiconductor  diode. 
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*2.6  Bipolar  Junction  Transistor  Model 

Weakly  nonlinear  effects  are  of  concern  primarily  when  the 
transistor  is  used  as  an  amplifier.  Consequently,  attention  in 
this  section  is  confined  to  applications  where  the  transistor 
operates  solely  in  the  amplification  region.  A  global  model  for 
the  bipolar  junction  transistor  in  the  amplification  region 
is  presented  in  Fig.  2.10.  The  terminals  denoted  by  E,  B ,  and  C 
are  referred  to  as  the  emitter,  base,  and  collector, 
respectively.  Loosely  speaking,  a  bipolar  junction  transistor- 
may  be  viewed  as  two  p-n  junctions  coupled  back-to-back.  In  the 
amplification  region  the  emitter-base  junction  is  forward  biased 
and  the  collector-base  junction  is  reverse  biased.  The  model  of 

f  \  ( 

Fig.  2.10  can  also  be  used  for  a  pnp  transistor  except  that  all 
variables  are  reversed  in  polarity.. 

The  seven  elements  in  the  global  model  are: 

1)  rb,  the  base  bulk  resistance 

21  Fje,  the  nonlinear  resistor  due  to  the  emitter-base 
junction.  Its  current-voltage  relationship  is  given  by 

lJE  ~  r8  <“*  <  >  -  «  •  ,2-6-1) 


( 
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Fig.  2.10  Global  model  for  the  npn  bipolar  junction  transistor  in 
the  amplification  region. 


3)  C.,_ ,  the  base-emitter  junction  diffusion  capacitance. 


4)  C  ,,  the  nonlinear  emitter-base  junction  diffusion 

U 

capacitance  whose  value  depends  linearly  on  the  emitter-base 
junction  current  in  accordance  with 


ca  --«4  **  •  <2-6'2’ 

5)  a ( i  ) Mi the  collector  current-controlled  current 
C  JE 

source.  Note  that  the  control  variable  i__  is  the  emitter-base 

JE 

junction  diode  current. 


6)  r  »  the  collector-base  junction  leakage  resistance. 


7)  C  ,  the  nonlinear  col lee tor -base  junction  space  charge 
capacitance  which  is  given  by 


C. 


K  V. 


rV 


CB 


(2.6-3) 


The  parameters  which  appear  in  the  controlled  source  and  in 
(2.6-1),  (2.6-2),  and  (2.6-3)  are  defined  as  follows: 


Xg  ■  reverse  saturation  current  of  the  emitter-base  diode, 

-19 

q  ■  electron  charge  ■  1.6  x  10  coulombs, 

-23 

k  ■  Boltzmann' 8  constant  «  1.38  x  10  joules/degree 
Kelvin, 

T  ■  junction  absolute  temperature  in  degrees  Kelvin, 
n  ■  emitter-base  diode  nonideality  factor, 

C  ■  emitter-base  junction  diffusion  capacitance  constant, 
o(Ic)  ■  common-base  normal-mode  dc  current  gain  ■  -1^/IE. 


'The  dependence  on  collector  current  is  given  by 

Xk 


1  ♦  +  a  ^lo  (  ) 

max 


(2.6-4) 


Cmax 


where  hpE  is  the  maximum  dc  current  gain  1 ^/I^,  Iq  is  the 
max 


B'  -max 


collector  current  at  which  *8  a  maximum,  and  a  is  a 

constant, 

M  ■  the  avalanche  multiplication  factor  given  by 


M 


[  i  -  <  )n  r1 

CBO 


(2.6-5) 
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where  is  the  avalanche  voltage  and  n  is  the  avalanche 
exponent, 

K  ■  collector-base  junction  space  charge  capacitance  at  1 
volt  collector-base  voltage, 
v  -  collector-base  junction  grading  constant. 


To  obtain  the  nonlinear  incremental  equivalent  circuit  for 
the  bipolar  junction  transistor,  it  is  necessary  to  expand  in 
power  series  the  incremental  currents  through  the  nonlinear 
elements  Rje,  Cd,  Cg,  and  a(ic)Mijg.  Since  NCAP  is  based  upon 
nodal  analysis,  all  of  these  expansions  must  be  in  terms  of 
incremental  voltages. 

The  incremental  current  through  may  be  written  as 

ije(t)  "  ^  9k  fv1e(t)]k  (2.6-6)  • 

where  the  linear  incremental  emitter  resistance  of  the  transistor 
is  defined  to  be 


_  1  nkT 

e  “g[  qT^T 


(2.6-7) 


and 


*1  Vl 


(2.6-8) 
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The  aeries  expansion  for  the  incremental  current  through  Cd  is 

lod<t)  “  Jcfo  K*T  "3t  {  Ivje(t>,kfl  }  (2.6-9) 

where 

Cd0  "  C2  (“IE)  (2.6-10) 

and 


Cdk  "  r  2  \  . 


(2.6-11) 


The  incremental  current  through  C8  is 

ics(t)  "  (c+T  ~3t  {  Jvd3(t)lk+1  }  (2.6-12) 

where 


c.o 

-  K 

v“^ 

CP 

(2.6-13) 

r 

-  (k  +, 

u  -  1) 

Lsk 

kVCB 

cs  (k-1)  ‘ 

(2.6-14) 

The 

final 

series 

to  be  presented  is 

that  for  the 
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current-controlled  current  source  a(icJMiJE.  This  is  obtained  by 
expressing  the  currents  ic  and  i^  in  terms  of  the  voltages  vCT 
and  vJE  .  It  follows  that  the  incremental  current  through  the 
controlled  source  may  be  written  as 

00 

i(t)  “  I  qk  tv.  (t)  ,  v.  <t)]  (2.6-15) 

8  k-1  K  ®  J® 

where  the  first  three  terms  are  given  by 

<Jl(Vdb'Vje>  '  “lIJEmlVdb  + 

^2(vdb'VJe’  '  <*lIJEm2  +  *2IJE^)vdb 

+  .  <Wl  +  .  2I0eWi*1>  v*vje 

+  +  *2nfri>  vj. 

«3(v*'V  ‘  +  2*2IJEmlm2  +  V^ 

*  <Wl  +  +  2V^1>  vrfjvje 

+  (.^gj  +  2“2Ij!f\)ml«2  +  2*2Vl«l)  Wj. 

♦  <”lV3  +  2*2IVl92  +  ‘ffab  vje  • 

(2.6-16) 


In  (2.6-16)  the  coefficients  o^,  c>2 ,  and  a3  are  associated  with 


(  th*  «(ic>  nonlinearity  of  (2.6-4).  They  are  given  by 


Ai.eo,  in  (2.6-16),  the  coefficients  mQ,  n^,  m2,  and  m3  are 
generated  from  the  avalanche  nonlinearity  of  (2.6-5).  They  are 


ml 

T 


When  the  transistor  is  biased  well  outside  of  the  avalanche 
region  (i.e.,  VCB«v^B0),  is  approximately  unity  while  m1#  m2, 
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4_  12  6-16)  involving 

.nd»3«.  — u*  Th*  * 

Vffl  can  then  be  cltclilt  foc  th.  npn 

The  nonli«.«  2.11.  ThU  *>d.l  i. 

;:r 


Fig.  2.11  Nonlinear  incremental  equivalent  circuit  for  npn 
bipolar  junction  transistor  in  the  amplification 
region. 


2.7  Junction  Field-Effect  Transistor  Model 


When  used  as  an  amplifier,  the  junction  field  effect 
transistor  is  biased  to  operate  in  a  region  of  its  static 
characteristics  known  as  the  saturation  region.  This  is  a  region 
where  weakly  nonlinear  effects  are  of  considerable  interest.  In 
the  saturation  region  a  global  model  for  the  n-channel  junction 
field-effect  transistor  is  shown  in  Fig.  2.12.  The  source,  gate, 
and  drain  terminals  are  denoted  by  s|  G,  and  D,  respectively. 


Fig.  2.12  Global  model  for  the  n-channel  junction  field-effect 
transistor  in  the  saturation  region. 

The  p-channel  junction  field-effect  transistor  has  the  same 
global  model  except  that  all  voltages  and  current  are  of  opposite 
polarity. 


The  four  elements  in  the  global  model  are 

1)  Rg,  the  aource  bulk  resistance. 

2)  i (v qj)  »  a  nonlinear  voltage-controlled  current  source. 
Note  that  the  control  variable  Vqj  is  the  voltage  be¬ 
tween  the  gate  terminal  and  the  internal  node  denoted  by 
<J.  Xn  the  saturation  region  an  analytical  expression 
•for  l(vGJ)  is 


2  +  4(^  - 


2fVCU+^7  , 


(2.7-1) 


3)  C„f  the  nonlinear  gate-to-source  capacitance  which 

Go 

is  given  by 


K  (-V  -  v  )' 
o  GT 


VGJ  <  0 


(2.7-2) 


4)  C(3D,  the  gate-to-drain  capacitance. 

The  parameters  which  appear  in  (2.7-1)  and  (2.7-2)  are  defined  as 
follows: 


IQ  ■  drain  current  at  maximum  dc  current  gain, 
max 

Vp  ■  pinchoff  voltage, 

♦  ■  JFET  barrier  potential, 

r  *  JFET  parameter  given  by 
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r  .  ^ 

Where  Ec  is  a  critical  field  value  aaaociated  with  the  nobility 
and  Lg  is  the  gate  length, 

V  ■  gate-to-source  capacitance  built-in  voltage, 

K  ■  gate-to-source  capacitance  for  “V0“vgj"1» 
m  ■  gate-to-source  capacitance  exponent. 


The  increnental  current  through  the  nonlinear  controlled 
source  i(VQj)  nay  be  expressed  as 


1(V 


I 

k-1 


*k 


vgj(t) 


k 


(2.7-4) 


To  sinplify  the  expressions  for  the  coefficients,  let 


1-e 


-r 


B 


A 


VGJ+* 

Vp  ^ 


(2.7-5) 
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junction  field-effect  transistor  is  shown  in  Fig.  2.13 


Fig.  2.13  Nonlinear'  incremental  equivalent  circuit  for  the  n- 

channel  junction  field-effect  transistor  in  the  satur¬ 
ation  region. 
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2.8  Vacuum-Tube  Diode  Model 


The  simplest  vacuum  tube  is  the  vacuum-tube  diode  which 
consists  of  two  electrodes  enclosed  in  a  vacuum.  A  global  model 
for  the  vacuum-tube  diode  is  presented  in  Fig.  2.14.  The  plate 


Fig.  2.14  Global  model  for  vacuum-tube  diode,  e  >0. 

B 

and  cathode  electrodes  are  denoted  by  P  and  K,  respectively.  The 
two  elements  in  the  global  model  arc 

1)  1^,  the  nonlinear  plate-to-cathode  resistance.  Its 
current-voltage  relationship  is  given  by 

3 

^  "  GeB  '  ®b  >  0  •  (2.8-1) 

2)  C  The  plate-to-cathode  capacitance. 

The  parameter  which  appears  in  (2.8-1)  is  defined  to  be 

G  ■  perveance. 

The  incremental  current  through  the  nonlinear  resistor  can 
be  written  as 


O 
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ir(t)  "  v 

r  k< 


where  the  linear 


z  \  [«b,t) 

“1  t 

incremental  resistance  of  the  diode  is 


(2.8-2) 


(.  ) 


3  -  J  x-i 


JL  »  (|G=S  > 

gl 


(2.8-3) 


and 


\ 


fl 


-1 


Vl 


(2.3-4) 


,  *.  „*reuit  for  the  vacuum-tube 

~  noniinear  incremental  eonivaU  t  .1  ^  ^ 

diode  is  presented  in  «9.  2'15’ 


circuit  consists  s 
resistor  rd. 


imply  of  the  capacitor  C* 
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in  parallel  with  the 


2.9  Vacuum-Tube  Triode  Model 


Assume  the  vacuum-tube  triode  is  biased  to  operate  'in  its 
amplification  region.  A  global  model  suitable  for  the 
amplification  region  is  shown  , in  Fig.  2.16.  The  grid,  plate,  an 
cathode 


Fig.  2.16  Global  model  for  vacuum-tube  triode  in  the  amplifi¬ 
cation  region. 


electrodes  are  denoted  by  G,  P,  and  K,  resDectively. 

The  four  circuit  elements  in  the  global  model  are: 

1)  Cgfc,  the  grid-to-cathode  capacitance. 

2)  C  ,  the  plate-to-gr id  capacitance. 

3)  C^,  ^e  plate-to-cathode  capacitance. 

4)  9(®g#«b)#  the  voltage-controlled  current  source 
representing  the  triode  space  current  which  is  given  by 
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•'E 


3 


(2.9-1) 


^(•a-V  *  Go  (1  ‘  r-2--1  <6g  +  ♦  +  F1  7 

<w 

The  parameters  which  appear  in  (2.9-1)  are  defined  as  follows: 

GQ  ■  perveance  for  eG"0, 

Ecinax  "  9ti^-to-cathode  voltage  for  perveance  equal  to  *ero, 

4>  *  offset  voltage  due  to  space  charge  effects, 
v  ■  triode  amplification  factor. 

Since  9(e_,e_  )  is  a  function  of  two  variables,  the 

G  B 
% 

incremental  space  current  is  obtained  from  a  two-dimensional 
Taylor  series  expansion  about  the  quiescent  operating  point.  The 
desired  expansion  may  be  written  as 


,(v^’  -  ^ v*  ‘w"2 

(n^  *  mj  /  0) 

The  e-pression  for  the  coefficients  is  simplified  by  defining 


G(eG>  - 

f  ^eG,eB^ 


(eG  +  <t> 


i. 

y 


3 

7 


(2.9-3) 


It  follows  that 
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_  q. 

J  l  \  (q^j)  l  »  1  r,(.a)  a2 

\>*2  "<«in  !m2M  “^Ty 

(qj  ♦  q2  -  »!> 


V.  "  Ec, 

S-'b  . 


(2.9-4) 


In  (2.9-4)  the  constraint,  g^+q  -m^,  indicates  that  only  terms 
for  which  the  indices  sum  to  m^  are  included  in  the  double 
summation.  The  first  several  coefficients  are  given  by 


*0,1 


Er  £*  i 

tl 

4nax 

t1- Ai  'eg  +  ♦  ♦  it1  7  -  et2- 

Snax  Snax 

tEG+*  +  ^)'7 

Snax 


»  3 

'eg  +  ♦  +  T1  7 


^  l1  '  «EG+*  +  V7 

Snax 


3G0  r_  ,  X  ,  Si.  7 

Wi -  [EG  +  ♦  +  V 

Snax 


'2,0  "  T 


-1  !*♦♦♦  V-sr2 


Snax 


Snax 


tEG  +  ♦  +  7 


(2.9-5) 


The  degree  of  a  term  in  (2.9-2)  is  the  sum  of  its  exponents. 
Henc*'  Kml>m2  is  the  coefficient  of  a  term  having  degree  (mi+ 
m2)  •  It  is  convenient  to  denote  by  gj<(eg»eb  )  the  sum  of  all 
terms  in  (2.9-2)  having  degrfee  k.  Then  (2.9-2)  may  be  expressed 


O 
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stVb> 


Z 

Jc— 1 


%.  '•a'V 


(2.9-6) 


where  the  first  three  terms  are  given  by 

°1(W  ■  K0,l  %  +  Kl,o  eg  5  %  *  *m  eg 

92<W  ‘  *0,2  ^  +  *1,1  «b  +  lt2,o  eg 

m  K0,3%  +  Kl,2*geb  +  k2,legeb  +  *3,0  ®g 

(2.9-7) 

The  nonlinear  incremental  equivalent  circuit  for  the  vacuum-tube 

triode  is  shown  in  Fiu.  2.17.  Note  that  KQn  is  the 

P 

reciprocal 


Fig.  2.17  Nonlinear  incremental  equivalent  circuit  for  the 
vacuum-tube  triode  in  the  amplification  region. 


of  the  triode  plate  resistance  while  K, =  g  is  the  triode 

10  ’n 
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2.10  Vacuum-Tube  Pentode  Model 


For  the  vacuum-tube  pentode,  attention  ia  also  focused  on 
the  amplification  region.  A  global  model  suitable  for  this 
region  is  shown  in  Fig.  2.18.  The  control  grid,  screen  grid, 
suppressor  grid,  plate,  and  cathode  terminals  are  denoted  by  , 


Fig.  2.18  Global  model  for  vacuum-tube  pentode  in  the  amplifi¬ 
cation  region. 

G2  ,  G3»  and  K,  respectively.  Since  the  suppressor  terminal  G^ 
s  typically  tied  to  the  cathode  K,  the  two  electrodes  are  shown 
as  a  single  terminal  in  the  figure.  The  interelectrode 
capacitances  are  represented  by  Cgk,  C^,  Cg^  ,  C  ,  ,  and 

Cpg  •  The  remaining  two  circuit  elements  in  the  global  model  are: 

1)  f  ® qi #  «Q2*  ®B)  *  the  voltage-controlled  current  source 
representing  the  space  current  between  the  cathode  and  screen 
which  is  given  by 
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fteG1'eG2'eB) 


•]  (er  +  4>  + 


G. 


3 

7 


1 

D 


i+  ( -^>m 

D  '  e„  1 


(2.10-1) 

2)  9  (®  gi'  *G2  '  ^  > '  the  voltage-controlled  current  source 

representing  the  space  current  between  the  cathode  and  plate 
which  is  given  by 


*<*G,'*G  'V  *  Go  [1  * 


■]  (e_  +  $  + 


3 

7 


<5^>m 

<=2 


max 


i+  th. 

n  e_ 

G, 


(2.10-2) 

The  parameters  which  appear  in  (2.10-1)  and  (2.10-2)  are  defined 
as  follows: 


G0  *  perveance  for  eQ1»0, 

Eaaax  "  contr°l  grid-to-cathode  voltage  for  perveance  equal  to 
zero, 

♦  ■  offset  voltage  due  to  space  charge  effects, 
u  •  control  grid  amplification  factor, 

D  ■  space  current  division  constant, 
m  ■  space  current  division  exponent. 

A  nonlinear  incremental  equivalent  circuit  for  the  pentode, 
which  is  appropriate  in  the  amplification  region,  is  obtained  by 
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ftxpsndin^  f  (®qi»  ^G2  *  #nd  9  »  ®G2*  ^g)  'In 

three-dimensional  Taylor  series  about  the  quiescent  operating 
point.  In  particular. 


f (e„l 


■,1-92-b'  ^  n2i„  n3l0  Anl-n2-n3  egJ  egl  eF  (2-111-3) 

nl«n2»n3  *  0 


and 


9(e„i /e^o/e.  )  * 

g!  g2  b  nl=0  n2»o  n3«0 


’  J  '  B„l'„2'„3  egJ  eg2  e?  <*-l*-<> 


nl*n2«n3  /  0 

where  the  coefficients  A  ir  ^ ,  n3  and  B^,  ^  ,  n3  are  given  by 


nl'n2'n3 


(nil)  (n2 1 )  (n3 f) 


a(nl.n2.n3,t(eGi>eG2teB) 

3eSi  8eSi  8eS3 


eGl*EGl 
®G2“EG2 
eB  “EB 


(2.10-5) 


and 

®dl*EGl 
®G2*EG2 
eB  "EB 

(2.10-6) 

Because  of  the  complexity  of  (2.10-5)  and  (2.10-6),  explicit 
expressions  for  the  coefficients  are  not  used  by  NCAP  for  their 


B 


nl'n2'n3 


3(nl+n2+n3)g(e 


(nl!J(n2!)(n3T T 


Gl'®G2'eB) 


9e 


nl  .  n2  -  n3 


3e":  9e 


G1  V~G2 


B 


evaluation.  Instead  the*  partial  derivatives  are  carried  out 
numerically.  Recall  that  Ani,  n2#  n3  and  Bnl,  n2  #  n3  are 
coefficients  of  terms  having  degree  (nl  +  n2  +  n3) .  It  is 
convenient  to  denote  by  fk(egl,  e  g2.  %)  and  gjc(^i  #  ^2  »  *b >  the 
sum  of  all' terms  in  (2.10-3)  and  (2.10-4),  respectively,  having 
degree  k.  Then  (2.10-3)  and  (2.10-4)  may  be  expressed  as 


f(*gl'*g2'*b)  "  fk(egl,eg2,eb) 


and 


£ 

k-1 


gk(egl,eg2,eb) 


(2.10-7) 


(2.10-8) 


tye  nonlinear  incremental  equivalent  circuit  for  the  vacuum-tube 
pentode  is  shown  in  Fig.  2.19. 
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Fig.  2.19  Nonlinear  incremental  equivalent  circuit  for  the 
vacuum-tube  pentode  in  the  amplification  region.  * 
Typically,  the- screen  grid  G  2  is  connected  externally  to  the 
cathode  through  an  impedance  which  behaves  as  a  short  circuit  at 
the  signal  frequencies  of  interest.  As  a  result,  e  G2  remains 
approximately  fixed  at  the  dc  value  E  q2  and  the  incremental 
voltage  eg2  is  approximately  zero.  Under  this  assumption,  the 
generator  f(eQ1  ,eG2  ,  eB)  is  incrementally  shorted  out  and  need 
not  be  considered  in  the  incremental  equivalent  circuit.  In 
addition,  the  generator  gfe^,  eG2  ,  63)  becomes 


J 


i 


G0  tl  -  E 


Br  Er  3 

X  <ea  **+.-%*  I 
Gi  u  1 


G2  • 

*  .  ,  o  ,m 

d+  (  e r } 


(2.10-9) 


The  incremental  current  corresponding  to  (2.10-9)  is  represented 
by  the.  two-dimensional  Taylor  series 
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(2.10-10) 

where  g^(e  )  denotes  the  sun  of  all  terms  in  (2.10-9)  having 

degree  k  and  the  coefficients  K  n^,  ^  are  given  by 


(n.+nj 

3  g(er  ,eJ) 


(2.10-11) 


Recognising  that  1)  C.  *  is  shorted  out,.  2)  C  is  in  parallel 
with  cg^g2  •  and  3)  Cg^  is  in  parallel  with  Cj*,  the  nonlinear 
incremental  equivalent  circuit  for  the  vacuum-tube  pentode 
becomes  that  shown  in  Fig.  2.20.  From  Fig.  2.20  it  is  apparent 
that 


1-75 


Pig.  2.20  Nonlinear,  incremental  equivalent  circuit  for  the 
vacuum-tube. pentode  in  the  amplification  region 
under  the  assumption  that  the  screen  grid  is 
incrementally  shorted  to  the  cathode. 

and  in  (2.2.10-10)  have  been  defined  such  that 

Sl^'V  *0,1  %  +  *1,0  agx  E  r“  %  +  egx  *  (2.10-12) 

Hence,  the  simplified  pentode  model  of  Pig.  2.20  is  identical  in 
form  to  the  triode  model  of  Fig.  2.17. 
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SECTION  3 

NONLINEAR  MODEL  .PARAMETERS  FOR  SEMICONDUCTOR  DEVICES  AND  VACUUM  TUBES 


3.1  Introduction 

In  order  to  use  NCAP  to  calculate  nonlinear  tranater 
functions  and  subsequently  determine  the  interference  to 
electronic  circuits,  a  detailed  circuit  model  must  be 
constructed.  Part  of  this  model  construction  includes  using  the 
stored  models  for  semiconductor  devices  and  vacuum  tubes  within 
the  computer  code  which  take  into  account  their  inherent 
nonlinear  behavior,  in  so  doing,  the  parameters  relevant  to  a 
given  device  model  must  be  determined.  This  section  will 
illustrate  how  the  device  parameters  can  be  obtained  so  that  a 
detailed  circuit  model  can  be  constructed. 

In  this  section,  nonlinear  model  parameters  for  a 

semiconductor  diode,  a  bipolar  junction  transistor  (BJT) ,  a 

junction  field  effect  transistor  (JPET) ,  and  vacuum  tube  diode, 

triode  and  pentode  will  be  determined.  The  techniques  employed 

in  determining  these  parameters  are  recommended  techniques,  i.e. 

they  are  not  the  sole  means  •  for  determining  the  required 

parameters.  However,  the  techniques  reported  upon  in  this 

section  have  been  used  with  a  sufficient  degree  of  confidence, 

12  3  4 

success  and  accuracy.  '  '  ' 

There  are  three  basic  techniques  for  obtaining  parameter 
values  for  the  models  used  in  NCAP.  The.  first  is  an  engineering 
guess  of  the  nominal  parameter  values.  This  could  yield 
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inaccurate  results  and  will  not  be  discussed  in  this  section. 
The  second  technique  employs  manufacturers  data  sheets.  The  use 
of  data  sheet  information  will  be  cited  where  applicable.  A 
major  drawback  with  this  technique  is,  in  general,  insufficient 
data  are  provided  to  determine  all  of  the  model  parameters.  The 
third  technique  is  the  measurement  technique.  This  will  be  the 
major  topic  of  discussion  in  this  section.  A  combination  of  dc, 
pulse  and  high  frequency  measurements  will  be  addressed  relevant 
to  the  parameter  values  for  the  various  models. 

It  is  appropriate  to  mention  at  this  point  that  there  are 
two  factors  relevant  to  parameter  value  determination  and, 
therefore,  circuit  model  construction,  which  are  not  explicitly 
accounted  for  in  the  stored  NCAP  device  models.  These  are  1) 
parasitic  effects,  such  as  bond  wire,  lead  inductances,  and 
fringing  capacitances,  and  2)  thermal  effects.  Parasitic 
components  can  be  added  externally  to  the  device  models.  Thermal 
effects  can  be  partially  accounted  for  by  modifying  appropriate 
parameter  values  to  reflect  actual  temperature  at  the  operating 
point. 


REFERENCES 


J.  w.  Graham  and  L.  Ehrman,  "Nonlinear  system  Modeling  anci 
Analysis  with  Application  to  Communication  Receivers," 
Technical  Report  RADC-TR-73-178,  Rome  Air  Development 
Center,  Griffiss  Air  Force  Base,  New  York,  June  1973.  (AD 
766  278) 

■J.  F.  Spina,  J.  Lexa,  D.  D.  Weiner,  "Computer  Modeling  or  a 
Solid-state  VHF  Tuner  Using  the  Nonlinear  Transfer  Function 
Approach,”  Technical  Memorandum,  Rome  Air  Development 
Center,  Griffiss  Air  Force  Base,  New  York,  1974. 

J.  C.  Bowers  and  S.  R.  Sedore,  "SCEPTRE:  A  Computer  Program 

«  . 

for  Circuit  and  System  Analysis,"  pp.  111-153,  New  Jersey, 
Prentice  Hall,  1971. 

J.  J.  Whalen,  C.  A.  Paludi,  Jr.  and  T.  F.  Fang, 
"Applications  of  the  Nonlinear  Circuit  Analysis  Program 
NCAP , "  1977  IEEE  International  Electromagnetic  Compatability 
Symposium  Record,  pp.  467-474,  Seattle,  Washington,  August 
1977.  (IEEE  Pub.  No.  770H  1231-0  EMC) 


1-79 


3.2  Nonlinear  Model  Parameter!  for  the  Semiconductor  Diode 

Figure  3.2-1  illustrates  the  model  for  the  semiconductor 
diode  used  by  NCAP.^  Additional  information  about  how  NCAP 
calculates  the  nonlinear  transfer  functions  using  this  model  is 
addressed  in  detail  in  Section  2.5. 


Fig.  3.2-1  Global  model  for  semiconductor  diode 

This  model  actually  represents  two  different  models;  one 
for  the  diode  in  the  forward  bias  region,  and  one  for  the  diode 
in  the  reverse  bias  region. 

The  parameters  required  by  NCAP  to  use  this  model  are  given 
in  Table  3.2-1.  Determination  of  these  parameters  will  be 
addressed  in  subsequent  subsections. 


TABLE  3.2-1 

NCAP  SEMICONDUCTOR  DIODE  MODEL  PARAMETERS 


NO.  NAME  UNITS 


DESCRIPTION 

dc  bias  current  diode 

diode  non-ideality  factor 
forward  biased  junction 
capacitance  extrapolated 
to  zero  current 
derivative  of  junction 
capacitance  with  re¬ 
spect  to  current 


V  magnitude  of  dc  reverse 

bias  voltage 

F  varactor  capacitance 

at  1  volt  reverse  bias 
varactor  capacitance 
exponent 

n  reverse  bias  leak- 

age  resistance 


a— -forward  bias 
b — reverse  bias 
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o 


3.2.1  Forward  Biased  Semiconductor  Diode 

A  typical  forward  biased  diode  it  illustrated  in  Fig.  3.2-2 

O - U — £>) - -O 

+  4 - V - ►  - 

Fig.  3. 2*- 2  Forward  biased  diode 


The  portion  of  the  global  model  given  in  Fig.  3.2-1,  used  by  NCAP 
in  the  forward  biased  region,  is  given  in  Fig.  3.2-3. 


Fig.  3.2-3  Global  model  for  semiconductor  diode  in 
forward  biased  region 


The  NCAP  parameter  ID,  dc  bias  current,  is  simply  the  dc 
current  at  the  operating  point  (0  pt.). 

The  diode  nonideality  factor,  n,  is  obtained  oy  measuring 
the  dc  forward  current-voltage  characteristic  which  is  governed 
by  the  mathematical  relationship  given  in  Eg.  3.2-1. 
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SVD 

r»  1  exp  (  5Et  >  "  1  1 


(3.2-1) 


In  the  forward  bias  region,  the  -1  term  can  be  neglected.  Baaed 
upon  Cq.  3.2-1,  a  semi-log  plot  of  Ip  vs  Vp  can  be  made.  Ip  is 
the  dc  forward  bias  current  plotted  on  the  ordinate  (log  scale) 
and  Vp  is  the  dc  forward  biased  voltage  plotted  on  the  abscissa 
(linear  scale).  A  typical  I-V  plot  is  illustrated  in  Pig. -3.2-4. 


Fig.  3.2-4  Typical  I-V  plot  for  determining  n 


The  simplest  means  of  determining  the  data  required  for  the 
plot  Ij,  vs  Vp  is  a  dc  measurement  shown  in  Pig.  3.2-5. 


Pig.  3.2-5  dc  measurement  set  up  for  determination  of  n 


A  transistor  curve  tracer  can  also  be  used  to  make  dc 
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measurements  to  obtain  the  appropriate  data  (see  Section  3.3.1). 

Since  the  diode  model  used  in  NCAP  does  not  explicitly 

account  for  thermal  effects,  measurement  techniques  for 

determining  parameter  values  should  be  chosen  such  that 

undesirable  thermal  effects . are  minimized.  When  dc  measurements 

are  used  to  determine  n,  erroneous  data  may  result  at  current 

levels  greater  than  10  mA.  This  is  illustrated  in  Pig.  3.2-4. 

To  overcome  these  effects,  a  pulse  measurement  scheme  may  be 
2 

used.  Figure  3.2-6  illustrates  a  pulse  measurement  scheme  used 

~  ,*ri 

to  obtain  the  necessary  data  (IF  vs  ^,)  to  determine  n  at  higher 
current  levels,  i.e.  greater  than  10  mA.  A  duty  cycle  of  10  to 
50  mS  (or  a  typical  pulse  duration  of  300  sec,  and  a  typical 
PRF  of  100  pulses/sec)  will  be  sufficient  to  obtain  the  required 


data. 


CURRENT 

PROBE 


Fig.  3.2-6  Pulse  measurement  set  up  for  determination  of  n 


It  should  be  mentioned  that  under  normal  (typical) 
conditions,  a  value  of  n  can  be  assumed  to  be  equal  to  1.0. 

The  diode  nonideality  factor,  n,  can  be  calculated  using  the 
data  plotted  in  Fig.  3.2-4  and  Eq.  3.2-2. 


1-84 


n 


(  q/kT  )  (  Vp2  -  Vpl  )•  /  In  (  Ip2/Ipl  ) 


(3.2-2) 


The  capacitive  effects  in  the  forward  biased  region  are  due 
to  the  nonlinear  diffusion  capacitance  of  the  p-n  junction,  cd . 
The  relationship  governing  this  effect  is  given  by  Eq.  3.2-3. 


(3.2-3) 


NCAP  requires  two  parameters  to. characterize  this  phenomena:  C ^ , 
the  forward  biased  junction  capacitance  extrapolated  to  zero 
current,. and  Cj,  the  derivative  of  junction  capacitance  with 
respect  to  current.  A  capacitance  bridge  can  be  used  to  obtain 
the  necessary  data  to  construct  a  linear  plot  of  C  vs  Ip.  Figure 
3.2-7  shows  the  diode  connection  on  a  Boonton  Capacitance  Bridge. 


Fig.  3.2-7  Diode  connection  on  capacitance  bridge 

A  plot  of  C  vs  IF  (on  linear  scales)  is  shown  in  Fig.  3.2-8. 
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Fig.  3.2-8  Typical  plot  of  C  vs  I 

The  value  for  C  ^  can  be  obtained  by  inspection  by  extrapolating 
the  data  to  ■  0.  cj  is  the  slope  of  the  straight  line  data. 


I 


(3.2-4) 


3.2.2  Reverse  Biased  Semiconductor  Diode 

A  typical  reverse  biased  diode  is  illustrated  in  Fig.  3.2-9. 

0 — - (<3 - o 

+  - - V - - 

Fig.  3.2-9  Reverse  biased  diode 


The  portion  of  the  global  model  given  in  Fig.  3.2-1,  used  by  NCAP 
in  the  reverse  biased  region,  is  given*  in  Fig.  3.2-10. 


o 


V0  <0 


Pig.  3.2-10  Global  model  for  semiconductor  diode  in  reverse 
biased  region. 


The  NCAP  parameter  VD,  magnitude  of  the  dc  reverse  bias 
voltage,  is  simply  the  dc  voltage  at  the  reverse  bias  operating 
point. 

The  reverse  bias  leakage  resistance,  R,  is  obtained  from  the 

dc  measurement  scheme  illustrated  in  Fig.  3.2-11. 

I00K 


AMMETER 


Pig.  3.2-11  Experimental  set-up  for  determining  reverse 
biased  dc  leakage  resistance  R 


A  characteristic  plot  (on  linear  coordinates)  can  be  made  from 
the  measured  data.  A  typical  plot  of  reverse  current  )  in  nA 
vs  reverse  voltage  (VR)  in  V  is  shown  in  Fig.  3.2-12. 
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*»  *> 


9  n  A 

|I„|  (A*) 


Fig.  3.2-12  Typical  plot  of  I R  vs  VR  for  determining  diode 
leakage  resistance,  R 

The  value  for  R  can  be  determined  by: 

AVR 

R  -  Tr-  (3.2-5)  . 

°  K 

The  NCAP  parameters  K,  varactor  capacitance  at  1  volt 
reverse  bias  and  y,  varactor  capacitance  exponent,  characterize 
the  nonlinear  transition  capacitance  in  the  reverse  bias  region 
as  indicated  by  Eq.  (3.2-6). 

Cs  “  K  (  ’VD  )_y  VD  <  0  (3.2-6) 

A  log-log  plot  of  capacitance,  C,  vs  reverse  voltage 
(magnitude)  ,  VR,  can  be  obtained  by  using  capacitance  bridge 
measurement  techniques  (see  Fig.  3.2-13). 
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Fig.  3.2-13 


Measurement  set  up  for  determining  reverse 
biased  diode  parameters  K  and  y 


Fig.  3.2-14  Typical  plot  of  C  vs  V  R 

The  curve  of  C  vs  ^  will  result  in  a  straight  line  when 
plotted  on  log-log  coordinates.  Using  Eq.  3.2-6: 


■  lo*10  C.  “  1O910  K  •  W  1O,10  1  VD  1 


(3.2-7) 


where  |vR|  from  Fig.  3.2-14  equals  I  VR|  .  The  value  K  can  be 
determined  by  inspection  of  Fig.  3.2-14.  The  value  y  the 
slope  of  the  C  vs  VR curve. 


1-89 


3 . 3  Nonlinear  Model  Parameters  for  the  Bipolar  Junction 


Tranaiator  (BJT) 

The  model  for  the  BJT  in  NCAP  is  referred  to  as  the 
nonlinear  T  model.  Figure  3.3-1  illustrates  the  linear 
incremental  T  model.  Prom  this  model,  the  nonlinear  behavior  of 
the  BJT  is  accounted  for  as  shown  in  Section  2. 


Fig.  3.3-1  Linear  incremental  T  model  for  BJT 

TABLE  3.3-1  lists  the  parameters  required  by  NCAP  to  characterize 
the  BJT.  Also  given  are  the  order  of  input  and  the  units  of  each 
parameter. 
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TABLE  3.3-1 


NCAP  BJT  Model  Parameters 


Number 

Name 

Description 

Units 

1 

n 

Avalanche  exponent 

2 

VCB 

Collector-base  bias  voltage 

V 

3 

VCBO 

Avalanche  voltage 

V 

4 

U 

Collector  capacitance  exponent 

5 

xc  . 

Collector  bias  current 

A 

6 

c  max 

Collector  current  at  maximum 

dc  current  gain 

A 

7 

a 

hpE  nonlinearity  coefficient 

8 

hFE 

max 

Maximum  dc  current  gain 

P'V** 

9 

K 

Collector  capacitor  scale  fac¬ 
tor 

10 

n 

Diode  nonideality  factor 

11 

°je 

Base-emitter  junction  space 

charge  capacitance 

P 

12 

C2 

Derivative  of  base-emitter 

diffusion  capacitance 

P/A 

13 

rb 

Base  resistance 

n 

14 

rc 

Collector  resistance 

ft 

15 

C1 

Base-emitter  capacitance 

p 

16 

C3 

Base-collector  and  overlap 

capacitance 

p 
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It  should  be  mentioned  that  certain  measured  data  have  signs 
that  depend  upon  whether  the  transistor  is  npn  or  pnp.  Table 
3.3-2  indicates  the  parameters  of  concern,  their  signs  when 
measured  under  normal  conditions  in  the  amplification  region,  and 
the  sign  convention  used  by  NCAP.  Note  that  NCAP  requires 
positive  values  for  these  parameters  independent  of  the  type  of 
transistor. 


TABLE  3.3-2 

Sign  Conventions  for  NCAP  BJT  Parameters 


Parameter 

Sign (npn) 

Sign (pnp) 

Sign  NCAP 

V 

♦ 

+ 

CB 

V 

♦ 

• 

+ 

CBO 

I 

+ 

— 

+ 

c 

I 

♦ 

- 

+ 

max 

3.3.1  n  - 

Diode  Nonideality  Factor 

and  r  -  Base 

Resistance 


The  diode  nonideality  factor  (n)  will  be  considered  first. 

The  value  for  r.  can  be  determined  from  data  required  by  n. 
b 

Several  techniques  for  determining  n  will  be  illustrated.  These 
techniques  will  result  in  a  characteristic  plot  L'sed  on  the 
current/voltage  relationship  given  in  Eq.  (3.3-1)  for  the 
base-emitter  junction  '  biased  under  normal  conditions  in  the 
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amplification  region ,  i.a.  tha  emitter-base  junction  ia  forward 
biased  and  the  collector-base  junction  is  reverse  biased. 


r  /  qV’«*E  \  1  1 

-IE  "  XES  loXp  {  ""rficT  )  ’  1  1 


(3.3-1) 


L 


Based  on  Eq.  (3.3-1).  a  plot  of  dc  emitter  current  (-1  )  on  a  log 

E 

ordinate  versus  the  dc  base-emitter  voltage  (V  )  on  a  linear 

BE 

abscissa  can  be  made.  A  typical-  plot  is  illustrated  in  Pig. 
3.3-2. 


NT 


Pig.  3.3-2  Typical  plot  of  -Ig  vs  VBS 


O 


Figures  3.3-3  and  3.3-4  illustrate  dc  experimental  set  ups 
for  obtaining  the  characteristic  plot  of  I£  vs  VB£. 
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Pig.  3.3-3  Experimental  setup  used  in  determination  of 
diode  nonideality  factor 


current 

SOURCE 


Pig.  3.3-4  dc  Measurement  for  determination  of  n 


In  Pig.  (3.3-3),  V  is  the  base  bias  supply.  The  variable 

BB 

resistor  acts  as  a  voltage  divider  and  is  used  to  adjust  V  . 

BE 

The  collector  bias  is  controlled  by  the  variable  collector  supply 
on  the  curve  tracer  (internally).  This  may  vary  with  V  due  to 

BE 

a  voltage  drop  across  the  internal  resistance  of  the  curve 
tracer.  V  must  be  kept  constant,  at  the  operating  point  (0 

CE 

pt.)  voltage,,  by  adjusting  the  collector  supply.  The  collector 

current  (I  )  is  displayed  on  the  curve  tracer  by  a  dot.  Since 

C 
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_IE  '  JB  +  V  *n,,  *B  is  much  less  than  ,  then  -IE  ■  Ic- 
Therefore,  -Ig  vs  VBE  can  be  plotted  as  shown  in  Fig.  3.3-2. 
Figure  3.3-4  illustrates  a  simplified  technique  for  obtaining  -IE 
and  V  directly. 

BE 

It  is  observed,  using  the  test  setups  in  Fig.  3.3-3  and 

3.3-4,  that  the  dc  measurements  of  log(-I  )  vs  V _ deviate  from  a 

E  BE 

straight  line  because  of  increased  junction  temperatures,  high 
level  injection  effects  and  the  fact  that  the  base-emitter 
junction  is  actually  a  distributed  R-C  network. 2' 7 ' 8Typically ,  dc 


■  .1  mA.  For  -l_  >  .1  mA,  a 

L 

pulse  measurement  scheme  is  used  to  obtain  the  required  data. 
Figure  3.3-5  shows  the  pulse  measurement  setup  and  the  results 


measurement  are  valid  up  to 


are  illustrated  on  Fig.  3.3-2. 


CURRENT 
R  PROBE 


Fig.  3.3-5  Pulse  measurement  scheme  for  determination 


of  -I  vs 
E 


A  duty  cycle  of  1  to  5%  is  typically  sufficient  to  obtain  the 
necessary  data. 

From  the  results  plotted  in  Fig.  3.3-2,  n  can  be  determined 
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n  ■ 


(q/kT)  (VBE2  -  VBE2>  /  In  dB2/IE1) 


(3.3-2) 


The  base  resistance  is  a  difficult  paraaeter  to  measure 
because  the  internal  base  node,  xx  ♦  1,  shown  in  Pig.  3.3-1,  is 
an  inaccessible  node.  The  simplest  way  to  obtain  r,  is  to  use 

D 

the  dc  base-emitter  characteristics  as  shown  in  Pig.  3.3-2.  The 
value  of  r^  simply  equals: 

rb  =  VQ/IEQ  (3.3-3) 

where  I  is  the  emitter  bias  current  and  V  is  obtained  from  the 
EQ  Q 

plot  of  -Ig  vs  V0E  (Fig.  3.3-2).  A  more  precise  value  for  rb  can 
be  determined  using  high  frequency  measurements,  i.e. 
s-parameters  and  y-parameters.  Using  y-parameters: 


11 


Vvin 


out 


(3.3-4) 


1 

*  rb  +  hfe  *e  (3.3-5) 

Solving  for  rb 


r 


b 


a t 


l 


hfe(nkT/q<"IE)> 


(3.3-6) 


where  h  and  y^  can  be  obtained  by  using  a  network  analyzer,  n 

is  the  diode  non-ideality  factor,  k  ■  Boltzmann's  Constant,  T  is 

temperature  (preset  at  300#K) ,  a  is  the  electronic  charge  and  -l_ 

is  the  dc  operating  point  emitter  current.  This  information  can 

9 

also  be  determined  from  s-parameter  data.  Manufactures  data 

sheets  often  provide  Smith  Chart  plots  of  s  11(si  ).  The  value 

of  r.  can  be  determined  from  that  value  of  the  plot  which  is 
o 

purely  resistive.  10 


3.3.2  I  -  dc  Collector  Bias  Current,  I 


-  Collector 


max 


Current  at  "Maximum  dc  Current  Gain,  a  -  F- 

FE 

Nonlinearity  Coefficient,  hp£  -  Maximum  dc 
Cur  rent*  T!a  in 


max 


The  parameters  Ic ,  Ic  ,  hpE ,  hpE 

max  max 

an  empirical  relationship  given  by  Eo.  3.3-7. 


and  are  related  by 


‘FE  M*X 


FE 


1  *  a  l°no  <  I~ 


CMAX 


(3.3-7) 


To  obtain  these  parameters,  a  characteristic  plot  must  be 
made  of  hpE  on  the  ordinate  versus  I q  on  a  log  abscissa.  A 
typical  plot  of  hpE  vs  Ic  is  illustrated  in  Fig.  3.3-6. 
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Fiq.  3.3-6  Typical  plot  of  hpE  vs  Ic 

There  are  three  basic  techniques  for  obtaining  the  necessary 

data  to  plot  hpE  vs  Ic<  The  first  technioue  employs  a  transistor 

curve  tracer.  The  common-emitter  output  characteristics  are 

displayed  on  the  curve  tracer.  Values  for  IB  and  Ip  can  be 

k  k 

determined  by  inspection  (see  Fig.  3.3-7).  Note  that  I„  's  can 

Bk 

be  determined  from  the  number  of  steps  displayed  on  the  curve 

tracer,  and  I  's  are  obtained  by  a  straight  line  extrapolation 
ck 

of  the  In  curve  to  the  I„  axis. 

3k  c 


WCI  1(1 


Fig.  3.3-7  Common-emitter  output  characteristics 

Since  h  ■  I  /L.  ,  h  vs  I  can  be  plotted  as  shown  in  Fig. 

FE  C 

3.3-6. 

A  more  direct  means  for  measuring  the  necessary  dc  currents 
is  illustrated  in  Fig.  3.3-8. 


Fig.  3.3-8  Experimental  setup  for  determination  of 
h  pE  nonlinear  parameters 

These  measurement  techniques  are  dc  techniques  and  are  valid 
for  low  current  levels,  i.e.  1^  <  1  mA.  For  high  current  levels, 
i.e.,  Ic  >  1mA,  a  pulse  measurement  technique  is  used.  This 
technique  is  illustrated  in  Fig.  3.3-9. 
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Pig.  3.3-9  Pulse  measurement  system  for  determination  of 


A  duty  cycle  of  1  to  5%  is  typically  sufficient  to  obtain  the 

t 

necessary  data. 

vs  I  as  shown  in  Pig. 

v 

3.3-6,  h  ,  and  I  can  be  determined  by  inspection;  the 

max  max 

value  for  (a)  has  to  be  calculated.  The  h  nonlinearity 

FE 

coefficient  (a)  can  be  calculated  in  two  ways,  first  at  the  Q 
pt.,  based  upon  the  value  of  I  ,  or,  over  a  given  range.  At  the  Q 

Pt.  ! 


From  the  characteristic  plot  of  h^ 


and  over  a  given  range,  say  (1^  ,  hpE1 )  and  (I^*  ^2*  where  XC1 
<  JC  <  TC2 
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^hFE2  "  hFEl^hFEl  loq10 *  I 


LC1 


CMAX 


.  .  2  ,  *C2  .  , 

)  “  hFE2  loglO(  I  ^ 


CMAX 


(3.3-9) 


3.3.3  n  ~  Avalanche  Exponent,  V  -  Collector-Base  Bias 

CB  __  _ 

Voltage,  V  -  Avalanche  Voltage,  r  -  Collector 

C  BU  _ C _  _ 

Resistance 

The  avalanche  effects  in  the  BJT  model  are  accounted  for  by 

an  avalanche  multiplication  factor  (M) .  The  value  for  M  is 

calculated  internally  by  NCAP.  However,  values  for  n#  V__  and 

CB 

V  must  be  provided.  These  parameters  can  be  determined  from 
CBO 

the  common-base  characteristics  of  the  BJT.  Figure  3.3-10 
illustrates  the  typical  common-base  characteristics  as  displayed 
on  a  transistor  curve  tracer. 

AVAL ANCHC 


Fig.  3.3-10  Typical  common-base  characteristics  of  BJT 

The  value  of  rc  can  be  determined  from  the  slope  of  the  lc 
curve  corresponding  to  the  dc  emitter  current  IE.  Denote  this 
slope  by  m.  Then 
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_1 

m 


(3.3-10) 


*c  " 

The  parameter  value  n  can  be  calculated  using  two  data 
points  (Icl  and  Ic2)  given  in  Pig.  3.3-10  and  Eg.  3.3-11. 


n 


log 


10 


Bi 


,  B1  -  rCK<°>  .  . 

B2  -  W01  ’ 


log 


10 


(  VCB1  ) 


CB2 


(3.3-11) 


where 


B, 


- 


CB1 


B 


2 


*C2 


r 


c 


The  value  for 
estimated  from  the  vs  VCB 


V  0  can  be  determined  two  ways. 


characteristics  as  shown 


3.3-10,  or  it  can  be  calculated  from  Eg  3.3-12. 


It  can  be 
in  Pig. 
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u 


CB 


CBO 


[  1  - 


W0)  i 

- 7 —  J 

lr  -  CB 
c 


(3.3-12) 


where  VC0  and  ^  are  values  measured  in  the  avalanche  region. 

3.3.4  u  -  Collector  Capacitance  Exponent,  K  -  Collector 
Capacitance  Scale  Factor 

Associated  with  the  collector-base  junction  is  a  varactor 
capacitor  effect  governed  by  Eg.  3.3-13. 


(  ) 


K  V, 


-M 


CB 


(3.3-13) 


Values  for  the  parameters  K  and  u  can  be  obtained  from  a 
log-log  plot  of  Cc  vs  |vCQ  |  using  Eg.  3.3-13. 


logi0  Cc  =  lo^io  K  -  "  ^lO^CB1  (3.3-14) 

Equation  3.3-14  results  in  a  straight  line  when  plotted  on 
log-log  coordinates.  K  is  the  value  of  capacitance  where  |vCB|  ■ 
1  volt  and  -u  is  the  slope  of  the  line.  Figure  3.3-11  shows  a 
typical  plot  of  Cc  vs  |  VCB  |  . 
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In  order  to  determine  K  and  \i ,  C  must  be  measured  as  a 

function  of  IV _ I.  when  an  RX  meter  is  used,  the  collector  and 

emitter  leads  are  connected  across  the  terminals  of  the  meter, 
and  the  base  and  emitter  leads  are  shorted  together.  The 
necessary  data  can  then  be  obtained.  When  a  capacitance  bridge 
is  used,  Fig.  3.3-12  illustrates  how  the  device  under  test  should 


be  connected  to  obtain  the  appropriate  data. 


Fig.  3.3-12  Device  connection  pattern  to  determine  NCAP 
parameters  K  and  n  on  capacitance  bridge 
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3.3.5  C.  -  Base-Emitter  Junction  Space  Charge  Capacitance 

J® _ _ 

C'  -  Derivative  of  Base-Emitter  Diffusion 
2 

Capacitance 

The  parameters  required  to  characterise  the  base-emitter 
junction  capacitance  effects  can  be  obtained  in  a  similar  fashion 
as  those  for  the  collector-base  junction.  Figure  3.3-13  shows 
the  connections  using  a  capacitance  bridge  to  obtain  a  plot  of  C^ 
vs-Iq.  The  mathematical  relationship  which  governs  these  effects 
is  given  in  Eq.  3.3-15. 


(3.3-15) 


where  C,  is  the  parallel  combination  of  C.  and  C^ .  When  c  vs 
i  je  d  2 

-Ig  ia  plotted  on  linear  coordinates,  the  plot  is  a  straight  line. 
Values  for  Cjfl  and  C’  can  be  determined  by  inspection.  CjG  is 
the  intercept  and  C£  is  the  slope  of  the  line  (see  Fig.  3.3-14). 


Fig.  3.3-13  Bridge  connection  of  DUT  for  determination  of 
base-emitter  capacitance  parameters 
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Fig.  3.3-14  Typical  plot  of  C2  vs  -I E 

It  is  necessary  to  mention  at  this  point  that  more  detailed, 
complex  and  in  general,  more  accurate  techniques  can  be  employed 
to  obtain  the  parameters  which  account  for  the  capacitive 
nonlinear  it ies  in  the  DJT  model  However,  those  techniques 
outlined  in  this  section  are  both  sufficient  and  accurate. 

3.3.6  -  Base-Emitter  Capacitance,  C3  -  Base-Collector 
and  Overlap  Capacitance 

Cx  and  C3  are  parasitic  capacitances  which  account  for  the 
leads,  bond  wires  and  package  of  the  BJT.  These  capacitances  are 
actually  external  to  the  model  for  the  BJT  and  play  a  significant 
role  at  high  frequencies  (i.e.  above  100  MHz). 
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There  is  no  direct  means  of  evaluating  each  parasitic 
component  separately.  Approximate  values  can  be  obtained  using 
techniques  addressed  by  Fang  7  or  by  adjusting  element  values  to 
obtain  a  best  fit  between  measured  and  calculated  data.  When  no 
data  are  available  to  determine  values  for  and  C 3#  and 

may  be  set  to  zero. 
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3 . 4  NCAP  Junction  Field-Effect  Translator  Model  Parameters 


The  RADC  Nonlinear  Circuit  Analysis  Program  NCAP  can  be  used 
to  calculate  nonlinear  effects  in  electronic  circuits  containing 
Junction  Field-Effect  Transistors  (JFET) .  In  this  section 
procedures  for  determining  the  NCAP  n-channel  JFET  (n-JFET)  model 
parameters  will  be  given.  Similar  procedures  can  be  used  for 
p-channel  JFETs.  As  a  specific  example  the  NCAP  JFET  model 
parameters  will  be  determined  for  a  2N3823  n-JFET  biased  at  Vffi  ■ 
-1.557  V,  VDS  «  +15  V,  and  ID  ■  40  UA.  We  first  review  the  NCAP 
JFET  analytic  model.  Then  we  describe  the  exper intents  performed 
to  determine  model  parameters.  We  also  discuss  the  data 
processing  techniques  used  to  obtain  numerical  values  for  the 
model  parameters. 

3.4.1  NCAPTotal  Equivalent  Circuit  for  a  JFET 

The  nonlinear  device  model  used  for  the  n-JFET  shown  in  Fig. 
3.4-1  is  referred  to  as  a  total  equivalent  model.1  3  The 
capacitors  Cl,  C2,  and  C3  shown  in  Fig.  3.4-1  are  constant 
extrinsic  linear  circuit  elements.  These  capacitors  are 
parasitic  elements  which  have  values  in  the  0.1  to  5  pF  range. 
An  intrinsic  device  model  for  the  JFET  is  shown  enclosed  within 
dotted  lines.  The  capacitor  CGD  is  the  gate-drain  capacitance 
which  is  modeled  as  a  constant  linear  circuit  element.  The 
resistors  .  Ps  and  RD  are  the  intrinsic  source  and  drain 
resistances  which  result  because  the  diffused  gate  region  does 


1-110 


not  completely  overlap  the  conducting  channel  and  because  the 
metallised  source  and  drain  contacts  are  not  perfect 
short-circuits.  These  resistors  are  modeled  as  constant  linear 
circuit  elements.  The  resistor  Rp  is  not  actually  included  in 

• 

the  NCAP  JPET  model  but  can  be  inserted  as  an  extrinsic  circuit 
element.  The  two  nonlinear  circuit  elements  are  the. 
gate-source  capacitance  CGS  (vGJ  )  and  the  dependent  current 
generator  i(vGJ),  where  vGJ  is  the  voltage  across  CGS  .  The  DC 
voltage  VGJ  is  related  to  the  actual  DC  terminal  voltage  VGS  by 


v 


GJ 


VGS  -  XdV 


(3.4-1) 


where  ID  is  the  DC  drain  current.  The  DC  gate  current  IG  is 
assumed  to  be  zero. 

The  computer  program  NCAP  contains  an  analytic  model  for  a 

JFET  which  is  well-suited  for  nonlinear  analysis.  At  DC  the 

dependent  current  generator  i(v  GJ  >  “*uals  1  DSAT<VGJ>  "h«re 

I DSAT^GJ^  *8  the  DC  dra*n  current  in  the  region  where  the  DC 

drain  current  is  independent  of  the  DC  drain  voltage  VDg.  The  DC 

saturated  drain  current  I^_.mis  given  by 

DSAT  J 


31. 


DSAT 


max 


t-p  +  -yp2  +  4  (  3  ■  n  +  j  n1,5  1 } » 


(3.4-2) 


(VGJ  ♦  T)  /  (Vp  +  V), 


(3.4-3) 
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p  ■ 


I 


(3*4-4) 


(1  -  exp(-D)  /  r2 

where  f  is  the  built-in  potential,  Vp  the  gate-source  cutoff 

voltage,  and  I_  and  r  parameters.  Equation  (3.4-2)  for  I 

MX  DSAT 

contains  four  parameters  (L.  ,r(or  p),  V,  and  V  )  which  can  be 

max  p 

determined  experimentally.  Note  that  the  voltages  VGJ,^  ,  and  Vp 
are  negative  numbers  for  the  n-JPET.  Also  note  that  either  the 
parameter  r  or  p  may  be  used;  however,  as  will  be  seen 
subsequently,  it  is  simpler  to  use  the  parameter  p. 

The  nonlinear  gate-source  capacitance  CGg  is  modeled  using 
the  expression 

CGS<VGJ>  *  K  I  -VGJ  -  Vo  ]-*  VGJ  <  o  <3.4-S> 

where  K  is  a  constant,  VQ  a  built-in  potential,  and  the  exponent 
m  a  constant.  Note  that  for  the  n-JFET  both  V~,  and  v  are 

GJ  O 

negative  real  numbers.  We  expect  the  built-in  potential  Vq  to  be 
in  the  range  -0.5  to  -0.7  V  for  a  silicon  n-JFET.  The  exponent  m 
is  0.5  for  an  ideal  abrupt  pn  junction  and  0.333  for  an  ideal 
linearly  graded  pn  junction.  We  anticipate  a  value  for  m  in  the 
range  0.3  to  0.5.  Equation  (3.4-5)  contains  three  parameters  K, 
V0 ,  and  m  which  can  be  determined  experimentally. 

We  may  summarize  the  NCAP  ana'vtic  model  parameters  in  Table 
3.4-1. 
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TABLE  3.4-1 


NCAP  JPET  Analytic  Model  Parameters 


Vqs 


c 

GD 

f,  V  ,r  orp, 
p 

K,  m,  VQ 
Rd.  Cl,  C2 .  C3 


I 


D 


max 


DC  gate-source  voltage 
source  resistance 
gate-dram  capacitance 
parameters  of  I (VG7 ) 
parameters  of  C^g  (Vqj  ) 
extrinsic  circuit  elements 
(not  intrinsic  model 
parameters) 


The  designation  of  V  as  the  pinch-off  voltage  in  Refs.  (2) 

n 

and  (3)  is  not  standard.  As  used  in  NCAP  V  is  equivalent  to  the 

p 

gate-source  cutoff  voltage  (cutoff)  where  (cutoff)  is  the 
value  o  f  VGS  that  makes  ID  «  0.  The  experimental  procedures  for 
measuring  the  model  parameters  will  be  given  in  the  following 
sub-sections . 


3.4.2  Experimental  Methods  for  Determining  the  JFET  Parameters, 


V  V  hi 


and  r  (or  p) 


How  NCAP  JFET  model  parameter  values  are  determined  will  be 
illustrated  for  the  2N3823  silicon  n-channel  JFET.  Both  DC  and 
Audio  Frequency  (AF)  measurement  techniques  are  used  to  measure 


the  five  JFET  parameters, 


Vo  ' 


and  r  (or  p  ) . 


We 


max 
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T 


assuire  the  built-in  potential  lies  in  the  range  -0.5  to  -0.7  V 

for  a  silicon  n-JPCT.  We  use  an  iterative  procedure  to  determine 

V  and  R  .  First  we  obtain  an  approximate  value  for  V 
n  S  d 

experimentally.  Next  we  determine  a  value  for  R  experimentally. 

Then  a  value  for  Vp  is  determined  using  another  experimental 

procedure.  We  also  determine  a  set  of  values  of  I  pgAj.  versus 

from  which  final  values  for  the  parameters  L.  ,  r  (or  P  )  ,  and  V 

max 

can  be  determined  using  a  curve-fitting  technique. 

A  DC  experimental  method  is  used  to  determine  a  value  for 

Fs  .  $  The  JFET  is  operated  in  the  linear  region  where  the  DC 

drain-source  voltage  Vpg  is  much  less  than  i  Vp  |  .  The  JFET 

eouivalent  circuit  shown  in  Fig.  3.4-1  may  also  be  used  for  DC 

analysis.  Fo  r  VDS  <<  |Vp|  we  replace  the  current  generator 

i (V  )  by  a  resistor  RCl  (\/, )  where 
GJ  GJ 

1/RQl  =  Gq  [1  -y'n  ],  (3.4-6) 

where  Cq  is  the  wide-open  channel  conductance.6  Again  we  mention 
that  the  parameter  Vp  is  the  gate-source  cutoff  voltage.  For 
Vpg  <<  |v  |  we  define  the  resistance 

"dso  "  -  rd  +  RS1  +  V 

Using  Eas.  (3.4-3),  (3.4-6),  and  (3.4-7)  we  obtain 

rdso  *  (rd  +  -V  *  G0  ~l  X>  ,3-4'81 

where  _ 

X  =  (1  -  -y  n  r1 

Ty  plotting  experimental  values  for  R^g  versus  values  for  X,  we 
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can  determine  the  sum  R f;  ♦  Rg  from  the  intercept  on  the  RjjS0 
axis. 

Shown  in  Pig.  3.4-2  is  the  experimental  system  used  to 
measure  the  resistance  RDS0  a  2N3823  silicon  n-JPET. 
Conventional  laboratory  DC  digital  voltmeters  and  ammeters  were 
used  to  measure  VDS,  VGS  and  ID.  Values  of  ID  were  measured  as  a 
function  of  VGS  at  VDS  -  0.1  V  and  at  VDg  ■  0.2  V,  and  values  for 
RDso  were  calculated  using  Eg.  (3.4-7).  These  values  are  given 
in  Table  3.4-2  .  The  average  value  for  RD50  is  also  given. 


TABLE  3.4-2 


Values  for  RDS0  Measured  at  VDg  <<|Vp|  for  a  2N3823  n-JFET 


VGS 

rdso  ^ 

Wn) 

RDSO*n) 

V 

at  VDg  -  0.1  V 

at  VDg  -  0.2  V 

average 

0.0000 

246 

247 

246 

-0.200 

281 

278 

279 

-0.400 

319 

317 

318 

-0.600 

370 

377 

374 

-0.800 

486 

494 

490 

-1.000 

658 

690 

674 

-1.200 

1176 

1234 

1205 

-1.400 

4762 

4348 

4555 

-1.600 

-  T - 

— r 
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AMMETER 


♦  DIGITAL 
VOLTMETER 


Pig.  3.4-2  Basic  test  circuit  for  measuring  the  n-JFET 
DC  characteristics 


(OHMS) 


O 


Fig.  3.4-3  Values  of  RDS0  versus  X  for  several  values  of 
Y  using  Vp  «  -1.70  V  and  VGJ  =  VGg 


1-118 


Shown  in  Fig.  3.4-3  are  experimental  values  for  R^  plotted 

versus  calculated  values  for  X.  To  calculate  values  for  X,  we 

use  Eq.  (3.4-3)  and  we  assume  that  V.  ■  which  is  equivalent 

GJ  CS 

to  assuming  that  the  I^Rg  voltage  drop  is  negligible  (See  Eq. 
(3.4-1)).  Note  that  for  < <  | vp|  »  the  *D  values  are  small, 
and  the  term  I  R,  is  also  small.  If  subsequent  calculations 

L  2d 

indicated  that  the  1^  value  was  not  negligible  compared  to  V  , 
we  would  use  an  iterative  calculation  procedure.  We  determined 
an  experimental  value  of  -1.7  V  for  by  measuring  the 

gate-source  voltage  that  made  less  than  1  yA.  Using 

V  ■  V _ ,  and  V  ■  -1.70  V,  we  calculated  values  for  X  for  three 

GJ  Q)  p 

values  of  the  built-in  potential  V.  For  values  of  X  <  4,  a 

straight  line  did  pass  through  several  data  points  for  each  value 
of  V  as  shown  in  Fig.  3.4-3  .  Also  given  are  the  intercepts  R^ 
♦  Rs  for  these  three  values  of  V  .  We  assumed  Rg  ■  RR  and  Rg  ■ 
(Rg  +  l^)/2.5  Values  for  Rg  are  also  given  in  Fig.  3.4-3  .  We 

observe  that  for  V  in  the  range  -0.5  to  -0.7  V  that  Rg  is  in  the 
range  32.5  to  42.5  fl.  We  selected  the  intermediate  value  Rg  ■ 
37  n. 

The  experimental  system  shown  in  Fig.  3.4-2  was  also  used 

to  measure  DC  values  of  Ips^-p  versus  .  These  values  are  given 

in  Table  3.4-3  .  The  values  given  for  I__-_  are  the  values  read 

DSAT 

on  the  ammeter  corrected  for  the  current  through  the  digital  VM 
used  to  measure  Vpg  .  Also  given  in  Table  3.4-3  are  values  for 
Vqj  calculated  using  Eq.  (3.4-1)  with  Rg  ■  37  n. 

The  system  shown  in  Fig.  3.4-4  was  used  to  make 

audio-frequency  (AF)  measurements  of  the  n-JFET  incremental 
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transconductance  gm  as  a  function  of  VGS.  By  plotting  Sr/Xdsat  va 
Vgj  ,  the  value  for  Vp  can  be  checked.  To  make  the  gm 
measurements  an  AP  generator  was  connected  as  shown  in  Pig. 

3.4-4  .  A  dual  channel  AP  voltmeter  was  connected  to  drain  and 
gate  nodes.  When  connected  in  this  manner,  the  gain  reading  on 
the  dual  channel  voltmeter  is  given  by  the  product  (-gm)  (1  kft). 
During  the  AF  gain  measurements,  the  digital  VM  used  to  measure 
VDS  was  disconnected.  The  AF  gain  measurements  were  made  at  an 
AF  frequency  of  1  kHz.  The  values  for  gm  are  also  given  in  Table 
3.4-3. 


TABLE  3.4-3 

Values  for  IDSAT,  VGJ,  and  gm  as  a  Function  of  VGS  at  VDS  -  15  V 


and  with  Rg  ■  37 

n 

VGS 

1 DSAT 

VGJ 

V 

mA 

V 

mmho 

0 

4.3 

-0.159 

3.47 

-.200 

3.55 

-0.331 

3.38 

-.400 

2.85 

-0.505 

3.16 

-.600 

2.10 

-0.678 

2.85 

-.800 

1.55 

-0.857 

2.57 

-1.000 

0.967 

-1.036 

2.16 

-1.200 

0.535 

-1.220 

1.67 

1-120 


-1.400 

0.170 

-1.437 

1.00 

-1.600 

0.019 

-1.600 

0.254 

-1.700 

0.003 

-1.700 

0.053 

pinch-off 

voltage  Vp  can 

be  estimated 

using  an 

approximate  JFET  relation 

'  -(Vp  '  VGJ>  t  n*  <3-4-9' 

where  n  is  a  parameter.  Using  the  values  for  1  Dsat  and  qm  9iven 
in  Table  3.4-3,  values  for  IDSAT  /9m  were  calculated  and  plotted 
versus  VGJ  in  Fig.  3.4-5  .  The  intercept  on  the  VGJ  axis  yields 
a  value  for  V  ■  -1.70  V.  This  value  for  the  gate-source  cutoff 

r 

voltage  Vp  is  in  good  agreement  with  a  value  determined  by 
measuring  the  DC  gate-source  voltage  VGS  that  reduces  IDSAT  t0 
values  less  than  1  \i A.  (2) 

Having  determined  values  for  Vp  and  Rs  and  assuming  ¥  is  in 

the  range  -0.5  to  -0.7  V,  we  calculated  values  for  the  other 

unknown  parameters  ID  and  p  (and  r  )  in  (3.4-2)  using  the 

max 

computer  program  VJFET  given  in  Appendix  3.4A.  Using  the  data 

given  in  Table  3.4-3  for  I  and  V  and  parameter  values  V  ■ 

DSAT  GS  P 

-1.7  V,  R  ■  37  n,  and  f  ■  -0.5  V  as  inputs,  the  computer  program 

b 

VJFET  yields  many  sets  of  values  for  the  parameters  ID  and  p 

max 

(or  r  )  as  outputs.  Using  several  of  these  computed  sets  of 

727  ;  '  " 

When  the  value  for  the  gate-source  cutoff  voltage  Vp  measured 
using  DC  techniaues  is  well  defined,  the  AF  measurement  procedure 
used  to  check  the  Vp  value  may  be  omitted. 
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DSAT/gm(V) 


parameter  values  for  In  and  p  (or  r  )»  values  for  I_0__  versus 

max  dsat 

VGS  were  calculated  using  Eqs.  (3.4-2)- (3 .4-4) .  The  interactive 
computer  program  WPET  given  in  Appendix  3.4B  was  used.  Shown  in 


Pig.  3.4-6  are  two  calculated  curves  of  I  DSAT  versus  V^g  which 
fit  the  experimental  values  quite  well.  Other  values  for  I  DSAT 


versus  VGS  calculated  with 
(or  r)  were  also  plotted,  but 


different  values  for  In  and  p 

max 


curves  through  these  calculated 


values  did  not  fit  the  experimental  values  quite  as  well  as 

curves  1  and  2.  We  observe  that  curve  1  with  I.  ■  13.95  mA 

umax 

and  p  *  0.530  (r  *  1.1303)  appears  to  fit  the  experimental  values 


better.  These  values  will  be  used  in  the  summary  of  parameter 
values  given  in  Table  3.4-6. 
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DSAT  (wA) 


Vft. (VOLTS) 
OS 


.  3.4-6  Experimental  and  calculated  values  for  IDSAT  vs  VGg. 

The  calculated  values  for  IDSAT  are  obtained  using 
Eqs.  (3.4-2)  and  (3.4-4)  with  V  *  -0.5  V,  V  -  -1.7  V, 


Rs  ■  37  n,  and  with  the  parameters  values  for  Ij 
and  p  (or  D  indicated 


max 
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3.4.3  Experimental  Methods  for  Determining  the  JFET 


Parameters  K,  m,  V,  and  CGp. 

3. 4. 3.1  Measurement  of  the  gate-drain  capacitance  CQD 

The  gate-drain  capacitance  CGD  can  be  determined  from  a 

2  6 

measured  value  of  the  JFET  common-source  y-parameter  yrs»  #  The 
y-parameter  y  can  be  determined  from  a  reverse  voltage  gain 

&  S 

measurement.  Such  a  measurement  can  be  made  using  the 
experimental  system  shown  in  Fig.  3.4-7.  The  parameter  y^g  at 
frequencies  near  1  MHz  may  be  expressed  as 


'rs 


-jfc>cr 


ss 


(3.4-10) 


and  the  parameter  C  .  is  given  by 

L  SS 


C 


rss 


C2  +  CGD' 


(3.4-11) 


Using  the  equivalent  circuit  given  in  Fig.  3.4-1  for  the 
JFET,  an  expression  for  C___  can  be  obtained  in  terms  of  the 

ITS  5 

voltage  gain  Av  of  the  auxiliary  amplifier  and  voltage  ratio 

V  /V.  measured  on  the  dual  channel  AC  voltmeter  shown  in  Fig. 
B  A 

3.4-7.  This  expression  is 


C 


rss 


3180  pF  „  1 

ttmhzI  ~ir 


x 


(3.4-12) 


where  the  signal  generator  frequency  f  is  in  MHz.  Equation 
(3.4-12)  is  obtained  by  noting  that  when  the  reactance  of  the 
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o 


Fig.  3.4-7  Experimental  system  used  to  determine  the  gate-drain 


capacitance  CGD.  The  n-JFET  reverse  voltage  gain  is 
measured  at  frequencies  near  1  MHz. 
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coupling  capacitor  C12  is  much  less  than  50  SI,  the  parallel 
capacitors  C  __  and  C2  (in  the  JPET  model)  and  the  50  si  resistor  at 

vjU 

the  auxiliary  amplifier  input  form  a  voltage  divider. 

The  measurement  procedure  was  as  follows.  The  DC  voltage 

applied  to  drain  was  ♦  15  V.  The  DC  voltage  V„_  applied  to  the 

gg 

voltage  divider  R1/R2  was  varied  and  measured.  The  ratio  R1/R2 
equaled  1.000  so  that  VGS  ■  vgg/2.  The  signal  input  to  the  drain 
node  was  adjusted  to  obtain  a  dual  channel  voltmeter  channel  B 
reading  in  the  correct  range.  The  auxiliary  amplifier  nominal 
voltage  gain  Ay  was  40  dB •  the  actual  voltage  gain  Ay  was 
measured  in  a  separate  experiment.  Measurements  were  made  at  0.5 
MHz  and  2.0  MHz.  Using  Eq.  (3.4-12)  and  the  measured  values  for 
Ay  and  VB/Vft,  values  for  Crsg  were  calculated.  These  values  are 
given  in  Table  3.4-4. 
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TABLE  3.4-4 


Values  for  C  versus  V_e  for  a  2N3823  JFET 

iSS  ub 

Measured  at  0.5  MHz  and  2.0  MHza 


VGS 

Creo  -  C2 
rss 

+  C  GD 

,V«/2> 

at  0.5  MHz 

it  2.0  MHz 

volts 

PF 

pF 

-0.016 

0.98 

0.94 

-0.500 

0.90 

0.88 

-1.000 

0.88 

0.84 

-1.500 

0.87 

0.82 

-2.000  b 

0.88 

0.82 

-3.000  k 

0.87 

0.82 

(3) 

We  cannot  separate 

C2  from  C  GD  using  this 

measurement 

technique.  We  set  C2  ■ 

0  and  CGD  -  Crss 

at  the  OP 

value  of  VRS. 

In  the  example  used  in  this  section,  the 

QP  values 

VGS  " 

-1.577  V,  V  ■  +15  V, 
DS 

and  I  •  40  u A . 

D 

The  larger 

value  of  CRD 

at  the  OP  is  0.87  pF.  The  value  C  GD  ■ 

0.87  pF 

includes  the 

actual  value  of  C__  plus  parasitic  capacitances  between  the  gate 
lead  and  drain  lead  due  to  the  socket  and  wiring.  This  value  is 
given  in  the  summary  of  2N3823  JFET  parameter  values  in  Table 
3.4-6. 


(3) 

a  The  JFET  is  being  operated  with  V  ■  +  15  V. 

DS 

b  The  JFET  is  cutoff  with  V  <  -1.70  V. 

GS 
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An  alternate  method  for  measuring  the  gate-drain  capacitance 

C _ is  illustrated  in  fig.  3.4-8.  Although  this  method  was  not 

GD 

used  for  the  2N3823  n-JFET  being  used  as  an  example  in  this 

section,  it  has  been  used  successfully  with  other  2N3823  n-JFETs . 

The  measurement  method  illustrated  in  Fig.  3.4-8  provides  a 

direct  measurement  of  the  capacitance  C  on  the  1  MHz 

GD 

capacitance  bridge.  This  is  a  significant  advantage. 

3.4.3. 2  Measurement  of  the  gate-source  capacitance 
parameters  K,  VQ,  and  m 

The  gate-source  capacitance  C  can  be  determined  from  a 

GS 

measured  value  of  the  JFET  common-source  y-parameter  y,  .  The 

is 

y-parameter  y^g  can  be  determined  from  an  input  admittance 
measurement.  The  parameter  yig  at  frequencies  near  1  MHz  may  be 
expressed  as 

^is  “  Giss  +  -^iss'  (3.4-13) 


and  the  parameter  C,  is  given  by 

iss 


Ciss  *  C1  +  C2  ♦  CGD  +  CQS, 


(3.4-14) 


where  Cl,  C2,  C  and  C  are  the  JFET  model  parameters  shown  in 

GD  GS 

Fig.  3.4-1. 

The  test  circuit  for  measuring  C.  is  shown  in  Fig.  3.4-9. 

iss 

A  0.047  yF  capacitor  was  connected  from  node  D  to  node  S.  This 
capacitor  acts  as  a  short-circuit  at  1  MHz  and  connects  the  drain 
and  source  of  the  JFET.  A  15  V  DC  power  supply  .provides  the 
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shown 


drain  bias.  The  drain  power  supply  DC  current  was  measured  with 
an  ammeter.  A  1  MHz  capacitance  bridge  was  used  to  make  input 
admittance  measurements.  Connections  were  made  from  the 
capacitance  bridge  to  the  gate  of  the  2N3823  JFET  and  to  the 
signal  ground.  The  capacitance  bridge  contains  an  internal  DC 
power  supply  which  was  used  to  provide  the  gate-source  bias 
voltage  VGS*  A  digital  voltmeter  was  connected  to  the  terminals 
on  the  capacitance  bridge  provided  for  monitoring  this  DC 
voltage.  Input  admittance  measurements  were  made  using  the 
balanced  mode  method  described  in  the  capacitance  bridge 
instruction  manual.  The  measured  values  are  given  in  Table 
3.4-5. 
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TABLE  3,4-5 


Manured  Value,  of  VGg.  13,  and  Ciss  and  Derived  Value,  of  VGJ, 
*D  *n',  Clas  "  CGD  ■  c2  for  lh«  2N3823  n-JFET. 


VGS 

13 

Ciss 

*D 

VGJ 

C  .  -  C 

iss  i 

Volts 

mA 

pP 

mA 

Volts 

pF 

-0.000 

2.96 

4.31 

2.95 

-0.109 

3.44 

-0.104 

2.66 

4.15 

2.65 

-0.202 

3.28 

-0.202 

2.36 

4.02 

2.35 

-0.289 

3.15 

-0.302 

2.16 

3.97 

2.15 

-0.382 

3.10 

-0.406 

1.86 

3.88 

1.85 

-0.474 

3.01 

-0.501 

1.66 

3.78 

1.65 

-0.562 

2.91 

-0.602 

1.46 

3.74 

1.45 

-0.656 

2.87 

-0.802 

1.06 

3.58 

1.05 

-0.841 

2.71 

-1.002 

0.66 

3.51 

.65 

-1.026 

2.64 

-1.250 

0.265 

3.37 

.250 

-1.259 

2.50 

-1.504 

0.065 

3.22 

.050 

-1.506 

2.35 

-1 . 754  a 

0 . 015b 

3.08 

.000 

-1.754 

2.21 

-2 . 005  a 

0.015b 

2.97 

.000 

-2.005 

2.10 

-2 . 517a 

0 . 015b 

2.83 

.000 

-2.517 

1.96 

-3 . 009a 

0 . 015b 

2.72 

.000 

-3.009 

1.85 

-4 . 002a 

0.015b 

2.60 

.000 

-4.002 

1.73 

-5 . 001a 

0 . 015b 

2.50 

.000 

-5.001 

1.63 

(4) 

(4) 

a  The 

2N3823  n- 

JFET  is 

cutoff 

for  V  < 
GS 

.5  mA. 

-1.70  V. 

b  The 

voltmeter 

current 

is  .01 
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The  data  presented  ir.  Table  3.4-5  can  be  used  to  obtain 
values  for  the  parameters  K ,  vQ,  and  m  in  the  Cq.  (Vqj)  expression 
which  is  Eq.  (3.4-5).  First  the  DC  voltage  VGJ  was  calculated 
using  Eq.  (3.4-1)  with  Rs  -  37  fi  .  The  values  for  ID  were 

calculated  by  subtracting  0.015  mA  from  the  values  for  the  drain 
DC  bias  supply  current  13  given  in  Table  3.4-5  since  13  equaled 
0.015  mA  when  the  2N3823  JFET  was  cutoff.  These  values  for  ID  do 
differ  slightly  from  the  values  for  IDgAT  given  in  Table  3.4-3, 
but  the  differences  do  not  significantly  affect  the  VGJ 
calculation.  Next,  we  subtract  the  known  value  of  C2  +  CGD  (0.87 
pF)  from  the  Ciss  values  to  obtain  values  for  Ciss  -  C2  -  CGD. 
The  values  for  ID>  VGJ,  and  Cigg  -  CGD  -  C2  derived  from  the 
measured  values  for  13,  V__  and  C.  are  also  given  in  Table 

3.4-5  for  a  2N3823  n-JFET. 

Using  Eqs.  (3.4-5)  and  (3.4-14),  we  obtain  the  expression 

Ciss  *  c2  -  CGD  *  cl  +  kIvgj  +  vorm  •  (3.4-15) 

which  now  must  be  fitted  to  the  values  of  C.  -  C2  -  C^  versus 

iss  GD 

v  .  We  know  the  value  for  V  lies  in  the  range  -0.5  to  -0.7  V 
GJ  o 

and  the  value  for  m  lies  in  the  range  0.3  to  0.5  for  a  silicon 

JFET.  We  use  as  initial  trial  values  V  »  -0.5  V  and  m  ■  0.5 

o 

(the  value  for  an  abrupt  pn  junction  gate).  Having  chosen 

initial  values  of  V  «  -0.5  V  and  m  ■  0.5,  the  next  step  is  to 

o 

select  two  sets  of  values  for  C  .  -  C2  -C„_  and  V  ,  (with  V__  > 

iss  GD  GJ  GJ 

V  ■  -1.7  V)  and  to  solve  Eq.  (3.4-15)  for  Cl  and  K.  Using  the 
P 

values  C  -  C2  -  C  -  3.44  pF  at  V  ■  -0.019  V.  and  C  -  C2 
iss  GD  GJ  iss 
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-  C  -  2.51  pF  at  V„,  -  -1.754  V  (the  closest  value  to  V_e  -V  ), 
GD  GJ  GS  p 

we  obtained  Cl  ■  0.80  pF  and  K  ■  2.13  pF  -  V  *5.  How  well  these 

initial  choices  for  V  ,  m,  Cl  and  K  fit  the  data  can  be 

o 

determined  by  plotting  the  values  of  CQS  ■  C^gg  -  C2  -  CGQ  -  Cl 

versus  values  of  IV  +  V  |  on  log-log  paper  and  drawing  a 

straight  line  of  slope  -  m  through  these  values.  Shown  in  Fig. 

3.4-10  is  such  a  plot  for  the  parameter  values  C2  +  C__  *  0.87 

pF ,  Cl  ■  0.80  pF  and  Vq  ■  -0.5  V.  A  straight  line  with  slope  ■ 

-0.50  fits  the  plotted  values  very  well  for  V  „  >  V  .  (Equation 

GS  p 

(3.4-15)  applies  when  the  JFET  is  not  cutoff.)  The  value  of  the 
parameter  K  for  the  straight  line  plotted  in  Fig.  3.4-10  is  2.20 
pF  -  v-5  which  is  within  5  percent  of  the  value  calculated 
previously.  We  could  use  either  value.  Since  the  initial 

choices  for  V  ,  m,  Cl  and  K  fit  the  C.  versus  V_,  data  quite 

O  1SS  bJ 

well,  no  additional  effort  was  made  to  fit  other  values  of  V  ,  m, 

o 

Cl,  and  K  to  the  Cigg  versus  VQJ  data.  The  parameter  values  VQ  ■ 
-0.5  V,  m  -  0.5,  Cl  ■  0.80  pF,  and  K  ■  2.20  pF-V*5  are  given  in 
the  summary  of  the  2N3823  JFET  parameter  values  in  Table  3.4-6. 
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SLOPE  *-0.50 


u 

© 

1.0 

U 

0.8 

1 

CM 

O 

0.6 

1 

0.4 

• 

• 

0.3 

cJ 

C2  -I-  C 


(pF) 

0.87 


*  experimental  values 


\j.c  v.h  u.b  1.0  2.0  4.0  6.0  10.0 

I  V0J  +  V0  | (VOLTS) 


Fig.  3.4-10  Values  of  CGg  versus  JvGJ  +  VQ | .  Values  for 

Ciss  and  VGJ  are  ta,cen  from  table  5.  Equation 
(3.4-1)  with  R  «  37  is  used  to  calculate 

S 

values  for  V__ 


TABLE  3.4-6 


O 


NCAP  Analytic  Modal  Parameter  Values  for  a  2N3823  n-JFET  Biased 

at  V  -  -1.577V,  V  -  +15  V,  and  I  »  40  yA. 

GS  Ub 


v 

.  _ 

1.557  V 

!D 

■ 

13.95  mA 

GS 

Bs 

m 

37  n 

max 

** 

i 

m 

1.1303  or  p 

CGD 

m 

0.87  pF 

% 

m 

-1.70  V 

rd 

m 

Rs  -  37  Q 

m 

-0.5  V 

Cl 

w 

0.8  pF 

K 

m 

2.20  pF  -  V 

Cl 

m 

0 

m 

at 

0.5 

C3 

m 

0 

vo 

* 

-0.5  V 
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Appendix  3. 4 A 


VJFET  Program  for  Determining  I  and  D  (or  D 

D  K 

max 

from  Measured  Values  of  I  and  V__. 

DSAT  GJ 


100  DIMENSION  CUR (9) , VOLT (9) 

110C  ORDER  PAIRS  OP  CURRENT  AND  VOLTAGE 

120  DATA  CUR/4 .30,3.60,2.83,2.10,1.52,0.96,0.50,0.177,0.018/ 
130  DATA  VOLT/. 00, -.20, -.40, -.60, -.80, -1.0, -1.2, -1.4, -1.6/ 
140C  INPUT  PARAMETERS 
150  PRINT* "VPI" 

160  READrVP 
170  PRINT?  "RHO" 

180  READ -RHO 
190  PRINT- "VBI" 

200  READ- VDI 

210C******************************************************* 

220  DO  110  J  -  1.8 

230  PRINT  1 

240  KST  *•  1  +  J 

250  DO  110  K  -  KST, 9 

260  AID1  -  CUR (J ) * . 001 

270  AID2  -  CUR (K) *. 001 

280  VGl  »  VOLT ( J )  -  CUR (J) *.037 

290  VG2  •  VOLT (K)  -CUR (K) *. 037 

300  PRINT  1,VG1,AID1,VG2,AID2 

310  VP  -  -1.7 

320  DO  100  I  -  1,100 

330  A1  ■  AM ( RHO , VGl , VP , VB I ) 

340  A2  •  AM  (RHO, VG2 , VP , VDI ) 

350  Bl  ■  DMDR (RHO, VGl , VP , VBI ) 

360  B2  -  DMDR (RHO , VG2 ,VP, VBI ) 

370  FPR  ■  - (A2*Bl  -Al*B2)/(A2*A2) 

380  ERROR  »  +AID1/AID2  -A1/A2 

390  IF(ABS(ERROR) .LT. 0.0001)  GO  TO  105 

400  EOFPR  -  ERROR/FPR 

410  IF ((RHO  -EOFPR) .LT. 0.0)  EOFPR-RHO/2. 

420  100  RHO  -  RHO  -EOFPR 

430  105  CUR1  ■  AIDl/3 ./RHO/Al 

440  CUR 2  -  AID2/3 . /RHO/A2 

450  110  PRINT  1, CUR1.CUR2, ERROR, FLOAT (I ), RHO 

460  1  FORMAT  (IX ,9G17. 4 ) 

470  STOP 
480  END 
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490C  ********************************* 
500  FUNCTION  AM (RHO, VG,VP,VBI) 

510  A  -  (VG  +VBI ) / (VP  ♦VBI) 

520  B  -  1./3,  -  A  ♦  2. *  (A**l . 5)/3 
530  AM  -  -RH0/2.  +SORT(RHO*RHO+4.*B)/2 
540  RETURN 
550  END 

560C********************************** 

570  FUNCTION  DMDR ( RHO , VG , VP , VB I ) 

580  A  ■  (VG  +VBI )  /  (VP  +VBI) 

590  B  »  1./3.-A  +2 . * (A**l . 5) /3 . 

600  DMDR  ■  -0.5  ♦  RHO/SORT(RHO*RHO  +  4 
610  RETURN 
620  END 


******** 

0*B)/2 . 0 


u 


Appendix  3.4B 


JFET  Program  for  Calculating  Values  of 

I  vs  V  . 

DSAT  GS 


100  DIMENSION  CUR  (10) ,VOLT(10) 

110  REAL  IDl , ID 2 , I DM AX , IDMXl , IDMX2 

120  CNST(X)  «  4.*((l./3.- (X+VBI ) / (VP+VBI ) )  +  (2 ,/3 . ) * ( (X+VBI ) / (VP+ 
1304  VBI ) ) **1 . 5) 

140  V  (Y)  «.5*(-RHO  +  (RHO*RHO  +  4  .*(  (1  ./3 (Y+VBI )/ (VP+VBI )  )  + 
1504  (2./3.)M  (Y+VBI)  /(VP+VBI))  **1.5))**.  5) 

160  9  PRINT • "NUMBER  OF  PAIRS" 

170  READ  1,NT 

180  PRINT- "CURRENT  VALUES  IN  MA" 

190  READ  1  CUR 

200  PRINT: "VOLTAGE  VALUES  IN  VOLTS" 

210  READ  1 , VOLT 

220  PRINT: "SOURCE  RESISTANCE  IN  K  OHMS" 

230  READ  1 ,  RES 

240  PRINT: "PINCH-OFF  VOLTAGE  IN  VOLTS" 

250  READ  1,VP 

260  PRINT: "BUILT-IN  VOLTAGE  IN  VOLTS" 

270  READ  l.VBI 
280  GAMMA  -  .5 
290  N  »  NN  -  1 

300**************************************************** 

310  DO  2  X  •  1,N 

320  PRINT  11 

330  IDl  ■  CUR { I ) *.001 

340  VG1  •  VOLT (I )  -  CUR (I ) *RES 

350  PRINT:"  VG1  IDSl" 

360  PRIN^  7 , I , VG1 , IDl 

370  PRINT  12 

380  PRINT  6 

390  KST  »  I  ♦  1 

400  DO  3  J  «  KST,NN 

410  ID2  -  CUR ( J ) * . 001 

420  VG2  ■  VOLT (J)  -  CUR(J)*RES 

430  RATIO  -  ID1/ID2 

440  X  -  2 . *RATIO* (RATIO  -  1.) 

450  Y  »  RATIO*RATIO*CNST (VG2 )  -  CNST(VGl) 

460  Al  ■  X* (2.*Y-X*CNST(VG2) ) 

470  A2  ■  Y*Y 

(J 


1-141 


480  IP ( (-A2/A1) .GE.0)  GOTO  5 
490  PRINT r "RHO  IS  IMAGINARY" 

500  GOTO  3 
510  5  CONTINUE 

52§«****** *********************************  ************ 

530  RHO  -  SORT (-A2/A1 ) 

540  IDMXl  «  IDl/ (3 . *RHO*V (VG1) ) 

550  IDMX2  ■  ID2/ ( 3 . *RHO*V (VG2)  ) 

560  DO  8  L  »  1,1000 
570  X  «  GAMMA 

580  GAMMA  -  SQRT( (1-EXP (-GAMMA) )/RHO) 

590  ERROR  »  10 . ** (AINT  (ALOG10 (ABS (GAMMA) ) ) -6 . ) 

600  IF(ABS(GAMMA-X) .LT. ERROR)  GOTO  10 
610  8  CONTINUE 

620  PRINT: "GAMMA  IS  UNDE RTERMINED  AFTER  1000  ITERATIONS" 

630**************************************************** 

640  10  PRINT  7, J,VG2,ID2, IDMXl, RHO, GAMMA 
650  3  CONTINUE 
660  2  CONTINUE 

670  PRINT:"  "PRINT* "  ": PRINT:"  ": PRINT:"  " 

680  PRINT: "AGAIN  (YES-1  NO-0)" 

690  READ  1,1 

700  IF(I.EO.l)  GOTO  9 

710  1  FORMAT (V) 

720  6  FORMAT  (9X , 3HVG2 , 8X , 4HIDS2 , 10X , 6HIDSMAX , 9X , 3HRHO, 10X , 5HGAMMA) 
730  7  FORMAT(1X,I2,3X,F8.4,3X,F10.7 , 3X ,E13 . 6 , 3X,F9 . 6 , 3X ,E13 . 6) 

740  11  FORMAT (IX, 70 (1H-) ) 

750  12  FORMAT (IX, 13 (2H  *) ) 

760  STOP: END 


J 


O 
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ILLUSTRATIONS 


Figure  3.4-1 
Figure  3.4-2 

Figure  3.4-3 

Figure  3.4-4 

Figure  3.4-5 

Figure  3.4-6 


Figure  3.4-7 


Figure  3.4-8 


Total  equivalent  circuit  model  for  the  n-JFET. 
Basic  test  circuit  for  measuring  the  n-JFET 
DC  characteristics. 

Values  of  R versus  X  for  several  values  of 
DSO 

4*  usinq  V  *  -1.7  V  and  V  ,  *  V  _. 

*  p  GJ  GS 

AF  test  system  used  to  measure  n-JFET  trans¬ 
conductance  g  . 

m 

Values  of  1  DSA'/^m  versus  Voj*  T^e  data 
given  in  Table  3.4-3  are  used  to  calculate 

values  for  I  /a  . 

DSAT/gm* 

Experimental  and  calculated  values  for  I  „ 

DSAT 

vs  V  .  The  calculated  values  for  I _ are 

GS  DSAT 

obtained  using  Eqs.  (3.4-2)  and  (3.4-4)  with 
4'  ■  -0.5  V,  Vp  -  -1.7  V,  Rs  -  37 


and  with  the  parameters  values  for  I  and 

D 

p  (orT)  indicated.  max 

Experimental  system  used  to  determine  the 
capacitance  C--..  The  n-JFET  reverse 
voltage  gain  is  measured  at  frequencies 
near  1  MHz. 

Test  circuit  for  measuring  C  _  .  The  1  MHz 

rss 

capacitance  bridge  can  be  biased  externally, 
and  the  external  bias  at  the  high  terminal 
can  be  monitored  with  an  external  voltmeter 
as  shown. 
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Figure  3.4-9 


The  1  MHz 


Test  circuit  for  measuring  Cigg. 

capacitance  bridge  contains  an  internal 

bias  supply  which  can  be  used  to  control 

V  .  The  V  value  is  read  on  an  exter- 
GS  GS 

nal  digital  voltmeter. 

Figure  3.4-10  Values  of  C  versus  IV  +  V  I.  Values  for 

GS  1  GJ  o' 

C  and  V  are  taken  from  Table  3.4-5. 
iss 

Equation  (3.4-1)  with  R  ■  3 7 jtj  is  used 

s 

to  calculate  values  for  V  . 

GJ 
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3.5  Vacuum  Tube  Diode 

There  are  still  many  communication  systems  in  use  that  use 
vacuum  tubes.  The  Nonlinear  Circuit  Analysis  Program  NCAP  can  be 
used  to  calculate  nonlinear  effects  in  electronic  circuits 
containing  vacuum  tube  diodes,  triodes,  and  pentodes.  In  this 
section  procedures  for  determining  the  NCAP  vacuum  tube  diode 
model  parameters  will  be  given. 

3.5.1  NCAP  Model  For _ the  Vacuum  Tube  Diode 

The  nonlinear  device  model  used  for  the  vacuum  tube  diode 

shown  in  Fig.  3.5-1  is  referred  to  as  a  total  equivalent  model. 

The  capacitor  Cpk  is  the  plate-to-cathode  capacitance  and  the 

resistor  R  is  the  nonlinear  plate-to-cathode  resistance.  It  is 
D 

assumed  that  the  vacuum  tube  diode  is  operated  at  values  of 


'B 


Fig.  3.5-1  Total  equivalent  circuit  for  the  vacuum  tube  diode. 

total  instantaneous  plate-cathode  voltage  cQ  low  enough  so  that 

the  total  instantaneous  plate  current  i  is  spaced-charge  limited 

B 


so  that 


o 


*  Go 


i.  5 
B 


0 


(3.5-1) 


where  G  is  the  perveance.  The  perveance  G  does  depend  upon  the 
cathode  heater  current  and  the  age  of  the  tube,  but  it  is  assumed 
to  be  independent  of  the  plate  voltage  efi  .  The  NCAP  vacuum  tube 
diode  parameters  are  summarized  in  Table  3.5-1. 

TABLE  3.5-1 

NCAP  Vacuum  Tube  Diode  Model  Parameters 


Parameter 


G 

E 


B 


C 


pk 


Description  &) 

Perveance 

DC  plate-cathode  voltage  at  operating 
point 

Plate-to-cathode  interelectrode 
capacitance 


3.5.2  Determining  the  Vacuum  Tube  Diode  Parameters  G  and_Efi 


The  vacuum  tube  diode  perveance  G  can  be  determined  from  the 
diode  dc  plate  current  I B  vs  dc  plate-cathode  voltage 
characteristics.  An  experimental  system  for  measuring  the  I  B  - 
E  characteristics  of  a  6H6  vacuum  tube  diode  (dual  diode)  is 

B 
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mm 


(j  shown 
in  Fig 


(  ) 


Fig. 


O 


I 


in  Fig.  3.5-2.  Typical  I  -  E  characteristics  are  shown 

B  B 

3.5-3. 


DC 

POWER 
S  UPPLY 


.5-2. 


Experimental  system  for  measuring  the  6H6  (dual 
vacuum  tube  diode)  I_  -  E_  characteristics. 
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\  ) 


Fig.' 3.5-3.  Typical  dc  plate  current  vs  plate-cathode  voltage 
characteristics  for  the  6H6  vacuum  tube  diode. 
(From  K.  R.  Spangenberg,  "Vacuum  Tubes,"  p.  5, 
McGraw  Hill,  New  York,  1948.) 
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Using  the  data  shown  in  Pig.  3.5-3  and  Eq.  (3.5-1),  values 

for  the  ratio  I  /E  ^ were  calculated.  These  values  are  given  in 
B  B 

Table  3.5-2.  It  is  seen  that  values  for  the  ratio  X  /E  and 

B  B 

hence  the  perveance  G  are  not  constant  but  decrease  as  the  dc 
plate  current  Ifi  increases.  Over  the  current  range 

TABLE  3.5-2 

Values  of  Perveance  G  vs  DC  Plate  Current  I„ 

O 

For  a  6H6  Vacuum  Tube  Diode 


IB  (A) 

eb  (V) 

G  ■  Ib(A)/Eb 
(A/V1*5] 

.004 

1.5 

.00217 

.008 

GO 

• 

CM 

.00171 

.012 

4.1 

.00145 

.016 

5.3 

.00131 

.020 

6.3 

.00126 

.024 

7.4 

.00119 

.028 

8.5 

.00113 

.032 

9.6 

.00106 

12  to  32  mA  the  perveance  G  has  a  value  .00126  A/V1*5  ♦  161. 

The  recommended  procedure  for  selecting  a  value  for  the 
perveance  G  is  to  determine  the  vacuum  tube  diode  dc  operating 
point  (I  ,  E  )  and  to  calculate  G  using  G  ■  I _  (A)  .*  En  (V)1*5  . 

B  B  B  B 
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I 

I 

I 

For  example  for  a  6H6  vacuum  tube  diode  having  a  dc  operating 

point  I.  •  24  mA  and  ED  ■  7.4  V,  the  appropriate  NCAP  vacuum  tube 
B  B 

diode  model  parameter  values  are: 

®B  *  7-«  V 
G  =  .00119  A/V1*5 


3.5.3  Determining  the  Vacuum  Tube  Diode  Parameter  CpJt 

An  examination  of  data  provided  by  manufacturers  of  vacuum 
tube  diodes  indicates  that  a  value  for  the  plate-to-cathode 
capacitance  CpJt  is  not  always  provided.  If  that  is  the  case,  the 
capacitance  Cp}c  can  be  measured  directly  on  a  1  MHz  capacitance 
bridge.  The  first  step  is  to  connect  the  appropriate  terminals 
of  the  vacuum  tube  socket  to  the  capacitance  bridqe  and  to  null 
the  capacitance  bridqe  or  to  record  the  capacitance  reading 
without  the  vacuum  tube  diode  inserted  in  the  socket.  Then  the 
vacuum  tube  diode  is  inserted  in  the  socket,  but  the  cathode 
heater  is  not  connected  (cold  cathode)  .  The  increase  in 
capacitance  yields  a  value  for  one  interelectrode  capacitance 
combination.  Dy  measuring  several  interelectrode  capacitance 
combinations,  it  is  usually  possible  to  estimate  the 
interelectrode  capacitance  between  two  specific  electrodes  such 
as  plate  and  cathode.  An  example  will  be  given  to  illustrate  how 
this  can  be  accomplished. 
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When  manufacturer ’s  data  are  available,  their  use  is 
recommended.  For  example  the  following  capacitance  data  are 
provided  in  the  RCA  Receiving  Tube  Manual  for  the  6AL5  vacuum 
tube  diode.  The  6AL5  twin  diode  is  widely  used  as  a  detector  in 
FM  and  television  circuits.  The  data  provided  are: 

Direct  Interelectrode  Capacitances: 

P*  No.  1  to  C*  No.  1,  Heater,  and  Internal  Shield 
P  No.  2  to  C  No.  2,  Heater,  and  Internal  Shield 

C  No.  1  to  P  No.  1,  Heater,  and  Internal  Shield 

C  No.  2  to  P  No.  2,  Heater,  and  Internal  Shield 

P  No.  1  to  P  No.  2 

*  P  -  Plate 
C  ■  Cathode 

A  simple  model  illustrating  the  major  interelectrode  capacitances 
for  the  6AL5  twin  diode  is  shown  in  Fig.  3.5-4. 


2 . 5x10  ~12  F 

2 . 5x10  ~12  F 
-12 

3.4x10  F 

-12 

3.4x10  F 

-12- 

.07x10  F 
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Q 


Fig.  3.5-4. 


Simple  intete 
tube  such  as 


lectrode  model  for 
the  6AL5. 


a 


dual  diode  vacuum 


Using  the  direct  interelectrode 
interelectrode  capacitance  values  g 


capacitance 
iven  in  Table 


values , 
3.5-3 


the 

were 


determined. 
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TABLE  3.5-3 


Values  For  the  Interelectrode  Capacitances 


C 

plkl 

C 

p2k2 

C 

Plp2 

C 

klhs 

C 

k2hs 


Plate  No.  1 
Plate  No.  2 
Plate  No.  1 
Cathode  No. 
Cathode  No. 


to  Cathode  No.  1 
to  Cathode  No.  2 
to  Plate  No.  2 

1  to  Heater-Shield 

2  to  Heater-Shield 


2.5x10  "12  F 
2.5x10  “12  F 
.07x10  “I2  F 
.9x10  ”12  F 
.9x10 “12  F 


The  values  for  the  plate-cathode  interelectrode  capacitances 

C  and  C  „  can  be  entered  as  part  of  the  set  of  NCAP 

plkl  p2k2 

parameter  values  for  the  vacuum  tube  diode.  The  other 
interelectrode  capacitances  must  be  specified  separately  as 
extrinsic  (passive)  components.  It  is  important  to  note  that 
other  parasitic  capacitances  such  as  these  associated  with  the 
tube  socket  and  leads  connected  to  the  socket  terminals  may  be 
important  at  frequencies  greater  than  100  MHz.  Important 
parasitic  capacitances  can  also  be  specified  separately  as 
extrinsic  (passive)  components. 


1-155 


3.6  Vacuum  Tube  Triode 


The  Nonlinear  Circuit  Analysis  Program  NCAP  can  be  used  to 
calculate  nonlinear  effects  in  electronic  circuits  containing 
vacuum  tube  triodes.  In  this  section  procedures  for  determining 
the  NCAP  vacuum  tube  triode  model  parameters  will  be  given. 

3.6.1  NCAP  Model  For  the  Vacuum  Tube  Tr iode 

The  nonlinear  device  model  used  for  the  vacuum  tube  triode 

shown  in  Fig.  3.6-1  may  be  called  a  total  equivalent  model.  The 

interelectrode  capacitors  C  C  ,  and  C  .  are  defined  in  Table 

gk  pg  pk 

3.6-1.  The  triode  space-charge  current  g(eG,eB)  depends  upon  the 
total  instantaneous  arid-cathode  voltage  e„  and  total 
instantaneous  plate-cathode  voltaae  e_.. 

O 
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Fig.  3.6-1.  Total  equivalent  circuit  for  the  vacuum  tube 
triode. 


The  expression  for  the  current  g(e_,  e_ )  used  in  NCAP  is 

G  B 


g(er,  eR)  =  G  [1  -  eG/Ec  ]  [eG  +  <fr  +  eR/u]  *  (3.6-1) 

max 


The  NCAP  model  parameters  G  ,  \i  ,  <t> ,  E  ,  E  and  E  are  defined 

o  C  G  B 

.  _  .  ,  .  ,  ,  max 

m  Table  3.6-1. 
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TABLE  3.6-1 


NCAP  Vacuum  Tube  Triode  Model  Parameters 


Parameter 
E  _ 


E 


B 


max 


gk 

pg 

pk 


Description 

DC  grid-cathode  voltage  at  operating  point 
DC  plate-cathode  voltage  at  operating  point 
Perveance  for  e„  ■  0 

b 

Grid-Cathode  voltage  for  zero  perveance 

Offset  voltage 

Triode  amplification  factor 

Grid  to  Cathode  interelectrode  capacitance 

Plate  to  Grid  interelecti ode  capacitance 

Plate  to  Cathode  interelectrode  capacitance 


3.6.2 


Determining  the  Vacuum  Tube  Triode  Parameters  E  ,  E  , 

- _  _ G _ B 

0  »  E  ,  <t>,  and  u 

_ _ umax 


Values  for  the  vacuum  tube  triode  model  parameters  E_  and  E_ 

are  determined  from  the  dc  operating  point  data  for  the  triode. 

The  model  parameters  G  ,  F_  .  <}>  ,  and  y  are  determined  from 

Lmax 

the  triode  dc  plate  current  I0  vs  dc  plate-cathode  voltage  Eg 
characteristics  as  a  function  of  the  dc  grid-cathode  voltage  EG. 
The  IB~CB  triode  characteristics  are  usually  provided  by  vacuum 
tube  manufacturers.  If  manufacturer's  data  are  not  available, 
the  T'i*EB  characteristics  can  be  measured  using  the  experimental 
system  shown  in  Fig.  3.6-2.  Shown  in  Fig.  3.6-3  are  IB  vs  E0 
characteristics  for  one  section  of  a  12AT7  dual  triode. 
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Fig.  3.6-3. 


DC  PLATE  VOLTAGE  EB(V) 


DC  I_  -  C_  characteristics  for 
B  B 

12AT7  triode. 


Prom  Eg.  (3.6-1)  the  dc  plate  current  I  -an  be  written  as 

B 


G(Eg)  IEg  +  <J>  +  EB/p] 


(3.6-2) 


G(EG>  “  Go  11  "  VEC  J 


(3.6-3) 


where  E  and  E  are  the  dc  grid-cathode  and  plate-cathode 

G  B 

voltages.  It  is  informative  to  plot  values  of  I  vs  E  as 

B  B 

shown  in  Fig.  3.6-4.  From  Ea.  (3.6-2)  the  expression  obtained 

for  l2/3  vs  E  for  constant  E  ,  y  and  $  is  that  of  a  straiqht 
B  B  G 

line. 


Ig/3  -  (G(EG)2/3  /y]  [yEG  +  y*  +  Eg] 


(3.6-4) 


For  constant  E^,  y  ,  and  4>  the  slope  S(EG)  and  Eg  axis  intercept 
1(E)  from  Ea.  (3.6-4)  are 


Slope  S(Eq)  =  G(EG)</J  /y 


(3.6-5) 


Intercept  I (EQ)  =  -(yEG  +  y$) 


(3.6-6) 


From 

Fig.  3.6-4 

it  is  seen  that  the  plots 

Of  I 

vs  E0  are 

straight  lines  for 

values  of  Ig^3 

>  .02  A 2/3  (I 

B  >  2'8 

mA)  . 

The 

slope 

S(E  )  of  the 
n 

straight  line 

segments  are 

given 

in 

Table 
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3’*  2'  If  tl1*  Straight  line 

v«  fB  tor  I  2/3  >  „2  2/3  e9”6ntS  thrOU9h  the  values  of  i  2/3 

for  the  ,  ,  .  are  ex*raPolated  to  i  ■ V3 fh  Ib 

eb  a*»a  intercept  r(P,  .  8  '  the  V»1  ,es 

°bU1"--  If  the  value  for  a  ls  9lVe"  T*bla  3.,-,  are 

Can  be  ^mated  tram  COnsta"t-  values  tor  „ 

U  =  41  <V  /  (-0  v 


(3.6-7) 


The  calculated  values  for  u  for  a  specific  range  of  E_  values 

corresponding  to  ae  are  also  given  in  Table  3.6-2.  As  the 

G 

grid-cathode  dc  bias  becomes  more  negative,  the  values  for  u 
initially  increase  and  then  decrease.  Since  the  computer  program 
NCAP  uses  a  constant  value  for  u,  a  single  value  for  u  must  be 
chosen.  Values  for  g*50,  40,  and  30  are  possible  choices. 

A  value  for  can  ^e  determined  using  the  following 

procedures.  First  Eo.  (3.6-5)  is  written  in  the  form 

G(Eg)  =  [u ’ slope  S ( EG ) ] 1 * 5  (3.6-8) 

Next  using  Eq.  (3.6-8)  and  the  values  for  the  slope  S(E  )  given 

G 

in  Table  3.6-2,  values  for  G(E  )  are  calculated  for  u  values  of 

G 

50,  40,  and  30.  (See  Table  3.6-2.)  The  resulting  values  of 
G(E_)  vs  E_  are  plotted  as  in  Fig.  3.6-5.  The  plots  are  straight 

G  G 

line  segments  for  E  <  -1  V  in  agreement  with  Ea.  (3.6-3).  By 

G 

extrapolating  the  straight  lines,  values  for  G  (G  -  axis 

o 

interceDt)  and  F.  _  (E  axis  intercept)  are  obtained.  These 

max  G 

values  are  listed  in  Table  3.6-3. 
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TABLE  3.6-3 


Values  of  E„  and  G-  vs  u  For  One  Section  of  a  12AT7  Triode 
Snax  ° 


M 

Ec 

max 

Go 

V 

mA/V 

50 

-7.9 

2.5 

40 

-7.8 

00 

• 

30 

-7.9 

1.15 
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PERVEANCE  G  (mA/V 


.t 


12  AT7 

JU.  50 

40 

30 

_  4 

tf) 

Sec. 2 

(J>  ( V  )  -.18 

-.225 

-.30 

> 

G^mA/V1*8)  2.5 

1.8 

1.15 

■s. 

< 

Er  (V)  -  8.0 

-7.8 

-8.0 

i  3 

- 

'■'MAX 

0 


0 


I - 1 _ I _ ^  ^ 

-2  -4  -6  -8 

DC  GRID -CATHODE  VOLTAGE  EG  (  V  ) 


Fig.  3.6-5. 


Perveance  G  vs  dc  grid-cathode  voltage  EG. 

the  G  -  axis  intercept*  F  is  the 

cmax 


intercept. 


Gq  is 


Eg  axis 


It  is  observed  that  the  value  for  Eg.  is  independent  of  the 

ul&X 

value  used  for  y .  However,  the  value  for  GQ  does  depend  upon  the 
y  value.  The  y  value  is  determined  by  the  triode  operating 
point.  For  example  if  EG  is  in  the  range  -1  to  >2  V,  the 
appropriate  value  is  y«*51.  (See  Table  3.6-2).  Then  values  for 
Gq  and  ♦  corresponding  to  y  ■  51  (or  50)  should  be  used  as  NCAP 
triode  model  parameters. 

The.  procedure  for  determining  NCAP  triode  model  parameters 
may  be  summarized  as  follows: 

(1)  Determine  EB  and  EG  from  the  triode  dc  operating  point 
data. 

(2)  Using  Table  3.6-2,  determine  a  single  value  for  y  based 
upon  the  value  for  E G  .  Interpolate  if  necessary. 

(3)  Using  Table  3.6-3  and  the  y  value  determined  in  step  (2), 

i 

determine  a  value  for  Ec  and  G  .  Interpolate  if 

max 

necessary. 

(4)  Use  Eq.  (3.6-6)  to  calculate  <t>.  To  obtain  a  value  for 
I(E^)  use  the  I(E~)  vs  E-  data  given  in  Table  3.6-2. 

li  \3  b 

Interpolate  if  necessary. 
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To  illustrate  the  procedure  a  specific  example  is  presented. 

The  12AT7  dc  operating  point  data  are  assumed  to  be  EB  »  +200  V 

and  Eg  ■  -1.5  V.  From  Table  3.6-2,  V  -  50.  From  Table  3.6-3 

E -  •  -7.9  V  and  G0  -  .0025  A/V1’5.  Using  Table  3.6-2,  the 

Sax 

intercept  I(EG)  is  estimated  to  be  70  V  at  EQ  ■  -1.5  V.  Using 
Eg.  (3.6-5)  with  I  -  70  V  and  EG  -  -1.5  V,  the  value  <fr  ■  +.1  V 
is  obtained.  (Note  that  the  value  obtained  is  positive.)  These 
$  values  are  summarized  in  Table  3.6-4. 

TABLE  3.6-4 

Values  of  NCAP  Model  Parameters  For  One  Section 
of  a  12 AT 7  Triode 


EB  »  +200  V 
»  -1.5  V 
U  ■  50 


-7.9  V 
.0025  A/V1,5 


<t>  ■  +.1  V 


3.6.3 


Determining  the  Vacuum  Tube  Triode  Parameters  Cg^., 

C  ,  and  C  . 
pg  px 


An  examination  of  data  provided  by  manufacturer's  of  vacuum 
tube  triodes  indicates  that  values  for  the  interelectrode 
capacitances  are  usually  available.  If  this  is  not  the  case,  the 
interelectrode  capacitances  can  be  measured  directly  on  a  1  MHz 
capacitance  bridge  using  procedures  similar  to  those  described  in 
Section  3.5.3  for  the  vacuum  tube  diode. 
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When  manufacturer's  data  are  available,  their  use  is 
recommended .  The  following  interelectrode  capacitance  data  for 
the  12AT7  twin  triode  are  taken  front  the  RCA  Receiving  Tube 
Manual . 


Direct  Interelectrode  Capacitance 


Without 

External  Shield 

Grid-Drive  Operation  pF 

Grid  to  Plate  (Each  Unit)  1.5 

Grid  to  Cathode  &  Heater  (Each  Unit)  2.2 


With 

External  Shield 
pF 
1.5 
2.2 


Plate  to  Cathode  and  Heater 


Unit 

Mo. 

1 

0.5 

Unit 

No. 

2 

0.4 

O 


1.5 


Pased  upon  the  direct  interelectrode  capacitance  data  provided, 

the  capacitance  values  given  in  Table  3.6-5  were  estimated.  The 

values  for  C  ,  C  ,  and  C  ,  can  be  entered  directly  as  NCAP 
gk  pg  pk 

parameter  values  for  the  triode.  If  a  shield  is  used,  the  plate 
to  shield  capacitance  CpS  can  be  specified  separately  as  an 
extrinsic  (passive)  component.  (The  shield  would  usually  be  at 
around  potential.)  Again  it  is  important  to  note  that  other 
parasitic  capacitances  such  as  those  associated  with  the  tube 
socket  and  leads  connected  to  the  tube  socket  terminals  may  be 
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important.  Important  parasitic  capacitances  can  be  specified 
separately  as  extrinsic  (passive)  components. 


TABLE  3.6-5 

Values  For  The  Interelectrode  Capacitances  for  the  12AT7  Triode 


C  ^  (Both  units) 
C  (Both  units) 

pg 

C  (Unit  1) 
pk 

C  ,  (Unit  2) 
pk 

C  (unit  1) 
ps 

C  (Unit  2) 
ps 


2.2  x  l<f12F 
-12 

1.5  x  10  F 
0.5  x  10"12F 
0.4  x  lP"1^ 
0.7  x  1P“12F 
1.1  x  10‘12F 
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3. 7_  Y??yyff,  Ty£®_??Dtode 

The  Nonlinear  Circuit  Analysis  Program  NCAP  can  be  used  to 
calculate  nonlinear  effects  in  electronic  circuits  containma 
vacuum  tube  pentodes.  A  schematic  for  the  pentode  is  shown  in 
Fio.  3.7-1.  A  pentode  is  normally  operated  with  the  suppressor 
and  connected  to  the  cathode,  and  the  NCAP  pentode  model  is 
restricted  to  this  case.  In  this  section  procedures  for 
determining  the  NCAP  vacuum  tube  pentode  model  parameters  will  be 
given. 


3.7.1  NCAP  Model  For  the  Vacuum  Tube  Pentode 

The  nonlinear  device  model  used  for  the  vacuum  tube  pentode 

in  Fiq.  3.7-2  may  be  called  a  total  equivalent  model.  The 

interelectrode  capacitors  C  ,  ,  C  .  C  *  C  ,  C  and  C  are 

'  gk  glg2  og  g2k  g2p  ok 

defined  in  Table  3.7-1.  The  pentode  space-charge  currents  f (c  , 

Vj  1 

e  ,  e  )  and  q(e  ,  e  ,  e  )  depend  upon  the  total  instantaneous 
G2  B  Gl  G2  B 

control-arid  cathode  voltaac  c_..,  total  instantaneous  screen-grid 

G  i 

cathode  voltage  c^_  ,  and  total  instantaneous  plate  cathode 

voltaqe  e  .  These  currents  arc  given  by 
B 


f(eGl®G2'eB) 


rax 


3 

Tl 


l 

D 


h 

G2 


(3.7-1) 
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Fin.  3.7-1.  Vacuuri  tube  pentode  schematic  illustrating  pentode 
total  instantaneous  terminal  currents  and  voltages 
(with  respect  to  cathode) .  In  normal  pentode 
operation  the  suppressor  grid  is  connected  to  the 
cathode  (e_,  *  P). 


The  MCAP  pentode  model  parameters  G  ,  E 

°  Cn 

E  ,  D,  and  m  are  defined  in  Table  3.7-1. 

B 


,  V  ,  *  .  eg1,  eG2> 


O 
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TABLE  3.7-1 


NCAP  Vacuum  Tube  Pentode  Model  Parameters  For  the  Case 
Where  Suppressor  Grid  is  Connected  to  Cathode 


Parameter 


Snax 


glg2 


Description 

DC  control-grid  cathode  voltage  at  operating  point 
DC  screen-grid  cathode  voltage  at  operating  point 
DC  plate  cathode  voltage  at  operating  point 
Perveance  for  eG1  *  0 

Control-grid  cathode  voltage  for  zero  perveance 
Offset  voltage 

Control-grid  amplification  factor 
Space-charge  current  division  constant 
Space-charge  current  division  exponent 
plate  to  control-grid  interelectrode  capacitance 
Control-arid  to  screen-grid  interelectrode 
capacitance 

Screen-grid  to  cathode  interelectrode  capacitance 

Sum  of  control-grid  to  cathode  interelectrode 

capacitance  C  plus  control-grid  to  suppressor- 

arid  interelectrode  capacitance  C  ?  C  »  C 

glg3  gk  glk 

+  C 

glg3 

Sum  of  plate  to  cathode  interelectrode  capacitance 
C  '  plus  plate  to  suppressor-grid  interelectrode 

PK 

capacitance  CpgJ.  C  p](.  Cpk  ♦  C^. 

Screen-arid  to  plate  interelectrode  capacitance 
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3.7.2  Determining  the  Vacuum  Tube  Pentode  Parameters  E 


r.f 


Values  for  the  vacuum  tube  pentode  model  parameters  E  , 

G  1 

F  ,  and  E  are  determined  from  the  dc  operating  point  data  for 
G2  B 

the  pentode. 


3.7.3  Determining  the  Vacuum  Tube  Pentode  Parameters  p  , 

E  r  ,  G  a  nd  <p  . 
rax  ° _ 

Py  connecting  the  screen-grid  G2  and  plate  P  together,  a 

pentode  can  be  operated  as  a  triode.  The  model  parameters 

parameters  G  ,  Ec  ,  <p  ,  and  D  can  be  determined  from  the  (Ip 
°  '“max 

I__  )  vs  F_  characteristics  as  a  function  of  the  dc  control-grid 
G2  B 

cathode  voltage  E^  where  I0  is  the  dc  plate  current,  lc2  the  dc 
screen-grid  current  and  E_  the  dc  plate  voltage.  Often  the 

O 

(I  +  i  )  vs  ED  characteristics  are  provided  by  vacuum  tube 

manufacturers.  If  manufacturer's  data  are  not  available,  the 

(I  4  I  )  vs  E_  characteristics  can  be  measured  with  the  same 
B  G2  o 

experimental  system  used  for  measuring  triode  characteristics 
shown  previously  in  Fig.  3.6-2.  The  triode  characteristics  for 
the  6AH6  pentode  are  shown  in  Fig.  3.7-3.  (The  plate  and 
screen-grid  are  connected  together  and  the  ammeter  A  reads 

xb  +  xg2  directly-) 

usina  Fgs.  (3.7-1)  and  (3.7-2)  the  sum,  of  the  dc  plate 

current  I„  and  dc  screen-grid  current  I__  can  be  expressed  as 

B  G4. 
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(IB  +  IG2)  "  G(EG1)  fEGl  +  ♦  +  E  r-/u] 


*Gl'  lbGl  +  ♦  +  EG2/Ul  (3.7-3) 


G(EG1)  "  Go  [  1  “  egi/ec  1 


(3.7-4) 
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o 


cathode  voltages.  The  next  step  is  to  plot  values  of 
(I  +  I__)  2>/3  vs  E__  as  shown  in  Fig.  3.7-4.  From  Eq.  (3.7-3) 

B  G2  G2 

the  expression  for  (Ifi  +  i  ^  2//3  obtained  is 

(IB  +  IG2>2/3  =  IG^GI1  2/3/ul  [ijEG1  +  ^  +  EG21  <3’7‘5> 


Intercept  I(EQl)  “  -(uEG1  +  U4> )  (3.7-7) 

From  Fig.  3.7-4  it  is  seen  that  the  plots  of  (I  ♦  I  )2^3  vs  E_- 

B  G2  G2 

are  straight  line  segments  in  agreement  with  Eo.  (3.7-5).  The 

slopes  S  (E  )  determined  from  the  straioht  line  segments  are 
Gl 

aiven  in  Table  3.7-2.  If  the  straight  line  segments  shown  in 

2/3 

Fig.  3.7-4  are  extrapolated  to  (IB  +  IG2)  ■  C,  the  values  for 
the  FG2  axis  intercept  I(EQ1)  given  in  Table  3.7-2  are  obtained. 
If  the  value  for  $  is  assummed  constant,  values  for  u  can  be 
calculated  from 

U  =  AI(EG1)/(-AEG1)  (3.7-8) 


The  values  for  the  ratio  y  calculated  using  Eo.  (3.7-8)  for  1  V 

changes  in  E  .  are  also  qiven  in  Table  3.7-2.  The  y  values  lie 
G1  M 
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in  the  range  26  to  34.  Since  the  computer  program  NCAP  uses  a 
constant  value  for  u  ,  a  single  value  for  m  must  be  chosen. 
Possible  choices  for  u  arc  35,  and  25. 
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A  value  for  Ep  can  be  determined  usina  the  followinq 

Hnnax 

procedure.  First  Pc.  (3.7-6)  is  written  as 


G(EG1> 


[  Vi  •  S(EG1)] 


(3.7-9) 


Text  using  Fo.  (3.7-9)  and  the  values  for  the  slope  S(Er,.)  given 


in  Table  3.7-2,  values  for  G(CG^)  are  calculated  for  u  values  of 

35,  3?  and  25.  (See  Table  3.7-2.)  The  resulting  values  of 

G(FGi)  vs  Eg1  are  plotted  in  Fig.  3.7-5.  The  plots  are  straight 

line  segments  in  agreement  with  Eo.  (3.7-4).  Py  extrapolating 

the  straight  line  searrents,  values  for  G  (the  G(E_.  )  axis 

o  Gl 

intercept)  and  *thc  CG1  a*i£  intercept)  are  determined. 

These  values  are  listed  in  Table  3.7-3. 


TAPLE  3.7-3 

Values  of  F  and  G  vs  u  For  a  6AH6  Pentode 
c  o 

max< 


M 

Fr  (V) 

'-max 

G  (irA/V 
o 

35 

-6.1 

5.1 

3  P 

-6.1 

4.1 

25 

-6.1 

3.1 

1-1  e3 


6 AH8  PENTODE 
TRIODE  -CONNECTED 


The  Er  value  is  observed  to  be  independent  of  the  value 
''max 

used  for  y  .  However,  the  value  for  GQ  does  depend  upon  the 
value  used  for  u  .  The  value  for  y  is  determined  by  the  pentode 
operating  point. 


The  following  procedure  for 
model  parameters  eg2>  Eg1<  Gq,  E 


determining  the  NCAP  pentode 

r  ,  y  ,  and  4>  is  recommended. 
Snax 


(1)  Determine  E G2  and  EG1  from  the  pentode  operating  point 
data. 


(2)  Using  Table  3.7-2,  determine  an  appropriate  value  for  y  . 
Interpolate  if  necessary. 

(3)  Using  Table  3.7-3  and  the  W  value  determined  in  step  (2), 

determine  the  values  for  Ec  and  G0-  Interpolate  if 

max 

necessary. 

(4)  Use  Eo.  (3.7-7)  to  calculate  4>  «  -I/y  -  EQ^.  To  obtain  a 
value  for  I(EG1)  use  the  I  ( E vs  EG1  data  given  in  Table 
3.7-2.  Interpolate  if  necessary. 


To  illustrate  the  procedure  a  specific  example  is  presented. 

The  6AH6  pentode  dc  operating  point  data  are  assumed  to  be  E  ,  * 

♦15P  V  and  E  *  -2  V.  Prom  Table  3.7-2,  an  appropriate  value 
G1 
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for  u  is  30.  From  Table  3.7-3,  the  values  for  E_  and  G  at 

Qnax  o 


30  are 


-6.1 


V  and  G  Q  ■  .0041  A/V 


1.5 


3.7-2,  the  intercept  I (En)  is  79  V  at  FQ1  »  -2  V, 


(3.7-7)  with 


Gl 


JG' 

-2 


V, 


30, 


From  Table 
Using  Eo. 
79  V,  the  value 


<t>  *  -I/p  -  E^  ■  -.33  V  is  obtained.  The  NCAP  pentode  model 
parameters  determined  in  this  subsection  are  summarized  in  Table 
3.7-4. 


TArLE  3.7-4 

Values  of  NCAP  Pentode  Model  Parameters  EG2,  EG^, 
u  ,  Fr  .  G_  ,  and  ^  for  a  6AH6  Pentode 


G2  * 

+  150  V 

Ec  «  -6.1  V 

iraxc  ^ 

Gl  * 

-2  V 

G  0  *  .0041  A/V 

W 

50 

*  =  -.33  V 

3.7.4  Determining  the  Vacuum  Tube  Pentode  Parameters  D 
and  m 


Py  operating  a  pentode  in  conventional  manner  with 
suppressor  grid  G3  connected  directly  to  the  cathode  K,  the 
pentode  parameters  P  and  m  can  be  determined.  Using  Eqs.  (3.7-1) 
and  (3.7-2)  the  ratio  of  the  dc  plate  current  IB  to  the  dc 
screen-grid  current  IG2  can  be  expressed  as 


IB/IG2 


■  D  [E./E-  n  ] 


m 


(3.7-10) 
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where  Eg  is  the  dc  plate  cathode  volta9e  and  EG2  is  the  dc 
screen-grid  cathode  voltage.  By  measuring  the  dc  currents  IB  and 
IG2  for  several  combinations  of  the  dc  voltages  Efi  and  EG2 

values  of  the  ratio  Ib/*G2  can  be  determined.  Shown  in  Fig. 

3.7- 6  is  a  logarithmic  plot  of  the  dc  current  ratio  I^/Ig2  vs  tbe 

dc  voltage  ratio  Ej/Eg2  for  .5  <  EB/EG2  <  3  for  several  values  of 

the  dc  control-grid  voltage  E  .  Ideally  the  logarithmic  plots 

G1 

of  I  /I  vs  E  /E  should  be  straight  lines  which  are 
B  G2  B  G2 

independent  of  E  but  a  weak  dependence  on  E  is  observed.  The 

G1  G1 

slope  of  the  log T I B  /Ir;2  1  vs  log(EB/EG2)  8trai9bt  line  segment 
yields  a  value  for  m  for  each  value  of  Efil .  The  value  of  IB  /I G2 
at  EB  /EG2  "  *  yields  a  value  for  D  for  each  value  of  EG1 .  The 

values  of  m  and  D  for  several  values  of  Eclare  given  in  Table 

3.7- 5. 


TABLE  3.7-5 

Values  of  m  and  D  for  a  6AH6  Pentode 


G1  <V> 

m 

D 

-.5 

.09 

4.0 

-1 

.11 

GO 

• 

n 

-2 

.13 

3.3 

-3 

.15 

3.0 
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rB /^62 


Fiq.  3.7-6.  Values  of  the  dc  current  ratio  I0/IG2  vs  eb/eG2  for 
several  values  of  the  control-grid  cathode  voltage 
Fg^  for  a  6/*H6  pentode.  The  suppressor-grid  is 
connected  to  the  cathode. 
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The  following  procedure  for  determining  the  NCAP  pentode 
model  parameters  m  and  P  is  recommended. 

(1)  Determine  the  appropriate  value  for  the  control-grid 
voltage  CG1  from  the  dc  operating  point  data. 

(2)  Use  Table  3.7-5  to  obtain  values  for  m  and  0.  Interpolate 
if  necessary. 

3.7.5  Determinino  the  Vacuum  Tube  Pentode  Interelcctrode 
Capacitances 


An  examination  of  data  provided  by  manufacturers  of  vacuum 
tube  pentodes  indicates  that  values  for  the  major  interelectrode 
capacitances  are  usually  available.  If  this  is  not  the  case  or 
if  additional  data  are  needed,  the  interelectrode  capacitances 
can  be  measured  directly  on  a  1  MHz  capacitance  bridge  using 
procedures  similar  to  these  described  in  Section  3.5.3  for  the 

vacuum  tube  diode. 

Interelectrode  capacitance  data  for  pentodes  are  usually 
provided  by  manufacturers  in  the  following  format. 
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Pircct  Interelectrode  Capacitances  (Approx.)* 


O 


Grid  No.  1  to  Plate  .02  pF 

Grid  No.  1  to  Cathode,  Heater,  Grid  No.  2  12  pF 

Grid  No.  3  and  Internal  Shield 
Plate  to  Cathode,  Heater,  Grid  No.  2,  4.7  pF 

Grid  No.  3  and  Internal  Shield 

*  Fxternal  shield  connected  to  cathode 


Usinq  the  direct  intereloctr ode  capacitance  data  given  above,  the 


values  C  j.  ■  .02  pf,  Cgk  «  12  pF  and  *  4.7  pF  are  obtained. 
The  values  for  Cgk,  C  ^  and  Cpg  can  be  entered  directly  as  NCAP 
parameter  values  for  the  pentode.  Other  interelectrode 
capacitance  values,  if  needed,  such  as  C  ,  Cg2k  and  Cg2p  would 
have  to  be  determined  experimentally  and  specified  separately  as 
extrinsic  (passive)  components. 


3.7.6  Summary  of  Pentode  Parameter  Values 


In  Table  3.7*6  the  NCAP  pentode  parameter  values  for  a  6AH6 
pentode  are  given  for  a  specific  dc  operating  point* 

Fg1  -  -2  V,  FG2  -  150  V,  and  Eg  -  300  V. 
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TABLE  3.7-6 


NCAP  Parameter  Values  For  a  6AH6  Pentode  At  A 


Specific  Operating 

Point 

«egi- 

v'  EG2  *  ;5"  V- 

and  E0  »  300  V) 

Parameter 

NO. 

Parameter  Nane 

Value 

units 

1 

Go 

.0041 

A/V  1 

2 

u 

50 

3 

D 

3.3 

4 

m 

.13 

5 

rGl 

-2 

V 

6 

Ecmax 

-6.1 

V 

7 

* 

-.33 

V 

P 

EG2 

150 

V 

9 

eb 

300 

V 

1C 

Cgk 

12x10  ~12 

F 

-12 

11 

cpg 

.02x10 

F 

12 

CDk 

4.7xlfl"12 

F 
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